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ABSTRACT 

We present an empirical dynamical model of the local interstellar medium 
based on 270 radial-velocity measurements for 157 sight lines toward nearby 
stars. Physical-parameter measurements (i.e., temperature, turbulent velocity, 
depletions) are available for 90 components, or one-third of the sample, enabling 
initial characterizations of the physical properties of LISM clouds. The model 
includes 15 warm clouds located within 15 pc of the Sun, each with a different 
velocity vector. We derive projected morphologies of all clouds and estimate the 
volume filling factor of warm partially ionized material in the LISM to be between 
~5.5% and 19%. Relative velocities of potentially interacting clouds are often 
supersonic, consistent with heating, turbulent, and metal-depletion properties. 
Cloud-cloud collisions may be responsible for the filamentary morphologies found 
in ~ 1/3 of LISM clouds, the distribution of clouds along the boundaries of the two 
nearest clouds (LIC and G), the detailed shape and heating of the Mic Cloud, 
the location of nearby radio scintillation screens, and the location of a LISM 
cold cloud. Contrary to previous claims, the Sun appears to be located in the 
transition zone between the LIC and G Clouds. 
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INTRODUCTION 



In their now classical theoreti cal mo dels for the int e rstellar medium (ISM). lField. Goldsmith. &: Habing 

fll969h . iMcKee fc Ostrikeil Jl977h . and IWolfire et all Jl995al jbl) assumed the ISM to be in 
thermal and steady-state equilibrium. In these models, three stable regimes co-exist in pres- 
sure equilibrium: the cold neutral medium (CNM) with temperature T > 50 K, the warm 
neutral (WNM) or ionized medium (WIM) with T ~ 8000 K, and the hot ionized medium 
(HIM) with T ~ 1, 000, 000 K. These models include heating by ultraviolet (UV) photons on 
grains and polycyclic aromatic hydrocarbon (PAH) molecules and cooling by various forbid- 
den lines and the hydrogen Lyman-a line. These models do not include gas flows or predict 
the expected sizes of the various components. Given the low density of interstellar gas and 
the presence of supernovae and strong stellar winds, one expects that the gas will b e far 
out of thermal and pressure equilibrium and be highly dynamic. Various reviews (e.g. 



Cox 



20051 : iMcCray fc Snowlll979l ) discuss these issues and highlight the complexity of the ISM. 



The Local Bubble (LB) is a region of low-density pre sumably hot gas exten ding in all 
directions to hydrogen column densities logiV(H I) = 19.3 (ILallement et al.ll2003l ). Its shape 
is determined by the onset of significant column density of Na I, indicative of a cold gas shell 
surrounding the LB. Although the LB is irregular in shape, it extends to roughly 100 pc 
from the Sun. For our purposes, we consider the local interstellar medium (LISM) to consist 
of the hot, warm, and cold gas located inside the LB. The LISM gas has been shaped by 
the supernovae explosions and winds of massive stars in the Scorpio-Centaurus Association 
and ionized and heated by radiat ion from hot stars and the Galactic UV background (e.g., 
Berghofer fc Breitschwerdtl |2002| ). and so should provide a useful test of interstellar gas 
properties in our Galaxy and the assumptions that underlie theoretical models of the ISM. We 
can now study the LISM in detail because the ground-level transitions of many neutral and 
ionized atoms present in the UV could be observed with the high-resolution spectrographs 
on the Hubble Space Telescope (HST). With resolutions as high a.s R = A/AA = 100,000 
(Af = 3.0 km s~^), both the Goddard High Resolution Spectrograph (GHRS) and the Space 
Telescope Imaging Spectrograph (STIS) have obtained stellar spectra containing numerous 
interstellar absorption lines. As described below, these ultraviolet spectra, together with 
ground-based spectra in the Ca II H and K resonance lines, provide critical data for sampling 
the kinematic and physical properties of warm interstellar gas along 157 lines of sight. 

The dynamical structure of the LISM has a direct infiuence on the structure of the 
heliosphere around our solar system and astrospheres surrounding other nearby stars. The 
extent of the heliosphere (astrosphere) is determined by the balance of momentum (pf^) 
between the outward moving solar (stellar) wind and the surrounding interstellar medium. 
Long-term variations in the solar wind strength are not well known, but observations of 
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astrospheres around young solar analogs provide clues as to what kind of wind the Sun had 
in its distant pas t. The solar wind 3 .5 billions years ago may have been ~35-fold stronger 
than it is today (IWood et al.ll2005al ). In contrast, density variations spanning 6 orders of 
magnitude are commonly observed throughout the general ISM. However, variations in the 
dynamical properties of the surrounding ISM can also cause significant variations in the 
structure of the heliosphere even between clo uds w ith little- t o-no d ensity variation. Reviews 



of heliospheric modeling include IZankI (119991 ) and 



Baranov 



1990f). and the det ection of as 



trospheres around nearby stars is reviewed by IWoodI (j2004j ). iMiiller et al.l (120061 ) explore the 
response of heliospheric models to various interstellar environments that exist in the LISM. 
Significant heliospheric (or astrospheric) compression can impact planetary albedos, atmo- 
spheric chemistry, and biol ogical mu t ation rates. Reviews of the implications of heliospheric 
variability are discussed by iRedfieldl (120061 ) and iFrischI (120061 ) . 



Crutcherl (119821 ) first noted that interstellar gas in the LISM flows in roughly the same 



direction away from the center of the Scorpio-Centaurus Association. iLallement fc Bertin 
(119921 ) then showed that the flow of interstellar gas in the direction away from the Galactic 
Center is consistent with a vector that differs somewhat from that of the gas in the Galactic 
Center direction. They coined the term AG Cloud for the former and G Cloud for the 
latter. The AG Cloud is now called the Local Interstellar Cloud (LIC), since the Sun has 
been presumed to be located just inside the LIC on the basis that the velocity of neutra l 
helium flowing into the heliosphere is consistent with the LIC flow vector (jWitte et al.lll993h . 
Lallement et al.l (119951 ) argued that the LISM has a comp lex velocity structure with at least 
seven clouds located within 12 pc. Using a larger data set, iFrisch. Grodnicki. fc Weltyl ((2002) 
were able to identify 7 clouds in the LISM on t he basis of their kinernatics. Using Na I spectra 
of stars in the Galactic anti-center hemisphere, iGenova fc BeckmanI (120031 ) identified 8 clouds 
of presumably cold gas lying beyond 50 pc of the Sun with velocity vectors very different 
from those identified in this paper. These clouds may lie at the edge of the Local Bubble or 
beyond. 

The present work expands on the earlier studies in two ways. First, we analyze a much 
larger data set consisting of 270 individual velocity components along 157 lines of sight 
through the LISM. Each velocity component provides kinematical information (i.e., radial 
velocity) of a parcel of gas that we can analyze together with other velocity components to 
identify velocity vectors and morphologies of different gas clouds. Second, high-resolution 
GHRS and STIS spectra for 55 of these velocity components allow us to measure the widths of 
absorption lines from at oms and ions of different a tomic weight to determine the temperature 
and turbulent velocity (IRedfield fc Linskyl|20g4b[). and for 6 5 of the velocity components D I 
(an excellent proxy for H I in the LISM; iLinsky et al.l 120061) is observed together wi th other 
ions which can be used to calculate metal depletions (cf., iRedfield fc Linskyll2004a] ). These 
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measurements allow us to determine some of the physical properties of the clouds in addition 
to their kinematical and morphological properties. 

In this paper, we use the term "cloud" to refer to a contiguous parcel of interstellar gas 
inside the LISM with homogeneous kinematical and physical properties. We determine the 
morphology of 15 such clouds (§ l2.3l and §[3]) by assuming that the interstellar gas flow inside 
each cloud is coherent and that the clouds have sharp edges. An upper limit to the distance of 
each cloud from the Sun is the distance to the nearest star whose spectrum shows a velocity 
component consistent with the cloud's velocity vector. We identify locations on the sky of 
possible cloud-cloud interactions and check whether these locations are consistent with other 
phenomena (§|5]). In subsequent papers we will describe how these cloud-cloud interactions 
could explain radio scintillation screens and the locations of cold clouds. The results may 
be used to test assumptions of sharp cloud boundaries, departures from coherent flow, and 
search for evidence of shear , cloud rotation and expansion, and alignment with magnetic 



fields fICox fc Heleniud 120031 ). 



2. DATA ANALYSIS 

2.1. Construction of LISM Observational Database 

The ability to assign large-scale dynamical flows to observed projected radial veloci- 
ties requires an extensive and densely sampled observational database. We have compiled 
the most comprehensive high-spectral-resolution observational database from absorption line 
transitions in the UV observed by HST and from the Ca II optical transition observed from 
the ground. Transitions in the UV are the most sensitive to the warm partially ionized clouds 
that populate t he LISM, while Ca II is the transition most sensitive to warm gas in the visible 



( iRedfieldl 120061 ). Only six sight lines have both UV and Ca II interstella r absorption detec 



tions : a Aql, a PsA, 6 Cas, rj UMa, a Gru, and e Gru (see references in iRedfield fc Linsky 



20021 ). Not all UV velocity components are detected in Ca II, but those that are agree in ob- 
served velocity very well. The absorption observed toward a Gru provides the only example 
where the ultra-high-resolution Ca II observ ations resolve multiple coraponents from ab- 



sorpt ion identified as a single UV component (ICrawford &: Dunkinlll995l : iRedfield &: Linsky 



20021 ). About 26% of the sight lines are observed in several ions. Multiple ion observations of 
the LISM along the same line of sight provide independent-projected velocity measurements 
and provide additional diagnostics of the physical properties of the material (e.g., depletion, 
temperature, ionization, etc). Even though only 1-3 velocity components are identified per 
sight line, moderately high spectral resolution is required to adequately separate and resolve 
individual absorbers with similar-projected velocities. 
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Our database, which includes 270 individual velocity components along 157 sight lines, 
is derived from: (1) the complete high-resol ution UV database of HST obser vations of LISM 
absorption toward stars within 100 pc (see iRedfield &: LinskyI 120021 . l2004al . and references 
therein), which represents 55% of the velocity components in our sample, (2) the h igh- 
resolution Ca II a bsorption measurements toward stars within 100 pc (IFrisch et al.ll2002l : and 
references within iRedfield fc LinskyI |2002| ). which represent 32% of the components in our 
sample, and (3) the moderate-r esolution UV database of HST observations of LISM absorp- 
tion toward stars within 100 pc (jWood. Alexander, fc Linskvll 19961 : IWood et al.ll2000l . l2005bl ). 
representing the remaining 13% of the components in our sample. Physical-parameter mea- 
surements (i.e., temperature, turbulent velocity, depletions) are available for 90 components, 
or one-third of the sample. 

All absorption is assumed to be caused by the LISM. Contamination of the absorption 
database by absorption caused by edge-on circumstellar disks, although possible, is highly 
unlikely. Not only are nearby stars with circumstellar material rare, the requirement of an 
edge-on orientation further limits the likelihood of observing such systems. Only a handful 
of such systems, which show circumstellar absorption, are known within 100 pc. The most 
prominent example is f3 Pic, whose spectrum shows stable absorption at the stellar rest 



frame and variable abs orption components, both due to circumstellar gas (iHobbs et al.lll985 



Brandeker et al.l 120041 ). in addition to a LISM component, resolved from the circumstellar 



material only in the he aviest (i.e., narrowest) ions, such as Fe II ( ILallement et al.l Il995 



Redfield &: LinskyI |2002| ). Onl y two other stars in the LI SM database have known edge-on 
circumstellar disks: (1) /5 Car (iLagrange-Henri et al.lll990l ) in which only the Na I absorption 
feature was observed to vary, whereas the Ca II absorption is relatively steady, and m atch 
the UV obse rvations jRedfield. Kessler-Silacci. fc Cieza 2007; Redfield fc Linsky 2002), and 
(2) AU Mic JKalas. Liu, fc Matthew? [20041 ) which shows no circumstellar abso rption in 



( Roberge et al.ll2005l : iFrance et al.ll2007l ) or other UV lines, including Lyman-a (IWood et al. 
2005bl ) and the single UV o bservation does not allow for any constraint on the constancy 
of the observed absorption (IRedfield &: LinskyI |2002| ) . Therefore, we have retained these 
absorption features in the LISM database, but their removal does not significantly change 
the velocity vectors determined for the clouds for which they are members. 

We focus on using Ca II because it primarily traces warm LISM gas, whereas Na I 



absorption toward nearby stars (e.g.. 


Blades, Wvnne- Jones. & Wavte 


1980 


1993; 


Weltv. Hobbs, & Kulkarni 


1994; 


Welsh et al. 


1994 


). Approximately c 



= 31%) of the sight lines in our sample also have Na I observations, of which LISM ab- 
sorption is detected along only 16 lines of sight. About one-half of the Na I detections are 
toward stars within 50 pc, so although relatively uncommon within 100 pc, Na I absorption 
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i s not significantly do minated by the cold gas located near the edge of the Local Bubble 



( iLallement et al.ll2003l ). Therefore, a Ca II absorption component that is associated with a 
Na I component at the same velocity is not necessarily indicative of distant gas. Practically 
all LISM Na I absorption components have companion absorption components in Ca II at the 
same velocity, which indicates that the cold gas detected by the Na I absorption is physically 
associated with the warm gas detected by the Ca II absorption, and not separate clouds at 
coincident velocities. It appears that much of the cold gas in the LISM is associated directly 
with warmer gas and these structures share a common velocity vector. 

Our combined database samples the sky unevenly because the sight lines were often 
selected to observe UV bright stars or for purposes other than measuring LISM absorption. 
Although this is the densest-sampled UV/optical database of LISM absorption to date, there 
remain significant regions of poor sampling in both space and distance. The distribution of 
sight lines is shown for all of our dynamical cloud structures in Figures [TT-fTSl The various 
symbols used to signify the sight lines, as well as a discussion of the morphology of the 
derived clouds, is provided in Section [31 The median angular distance from one sight line to 
its nearest neighbor is 6.6 degrees, ranging from observations of binary stars with angular 
separations of ~10 arcsecs (e.g., a Cen A and B, a CMa A and B), to the poorly sampled 
region near I = 137° and b = 49° where the maximum nearest neighbor separation is 21.5 
degrees. Areas of poor sampling limit our ability to detect dynamical cloud structures, as 
many sight lines through the same collection of gas are required to determine an accurate 
velocity vector. Poor sampling also limits our ability to estimate distances to structures in 
these regions. 



2.2. Criteria for Identifying an Interstellar Cloud 



We began our search for identifi able structures in th e LISM with th e properties of the 
Lie a s our prototype. As shown by ILallement fc Bertinl (119921 ) and by iRedfield &: Linsky 
(I2OOOI ). the Lie moves as if it were a rigid structure, that is the observed radial velocities 
toward nearby stars over a wide range of Galactic coordinates are consistent with a sin- 
gle velocity vector. The scatter of the measured radial velocities about the mean vector 
is generally l ess than 1 km s~^, which is similar to the absolute velocity precision of STIS 
echelle data. iRedfield &: LinskyI (120001 ) constructed a three-dimensional model for the LIC 
based on absorption-line data for 32 lines of sight. The edge of the cloud was determined by 
the measured H I column density along each line of sight and the assumption that the H I 
number density is the same throughout the LIC. This simple assumption cannot be readily 
tested and could be far from the truth. Thus the true shape of the LIC is not well deter- 
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mined and the question of its edge, whether it be sharp or gradual, is unknown. Although 



Redfield &: LinskyI (120001 ) concluded that the LIC is roughly spherical in shape, the shape of 



interstella r clouds often appea rs to be filamenta ry based on ari abun dance of observations 



Graham et aPllQQsh filamentary structur e 



of nearby (iFrisch fc Yorklll983l ) and distant (e.g. 
in the ISM, presumably organized by magnetic fields (I Jackson. Werner, fc Gautierl l2003l ). 
We have therefore not assumed any a priori shape for the clouds in the LISM. In practice, 
we have followed a few simple rules in identifying interstellar clouds. Figures [T]-[T5l show 
the spatial distribution and projected boundaries of the resulting dynamical clouds, and 
Tables [T]-[T5l list sight line membership of the clouds. 



We determine the three-dimensional heliocentric velocity vector (three free parameters: 
the velocity magnitude [Vq] and the direction in Galactic coordinates [Iq, 6o]) that best 
fits the radial velocity database, where K- = Vo(cos h cos 60 cos(/o — + &o sin h). Here 
Vo is >0 for downwind directions and <0 for upwind directions. The first application of 
this procedure to the entire database yields a velocity vector consistent with absorption 
due to the LIC. Since LIC absorption is seen over much of the sky, the LIC should 
have the greatest number of observed sight lines and dominate the dynamical fit of 
the entire database. We next delete the velocity component that most significantly 
disagrees with the predicted projected velocity for LIC absorption and then recompute 
the velocity vector that best fits the remaining points. This procedure is continued until 
a satisfactory fit to the data is derived. Our criterion for ending the iteration process 
is that the removal of next most discrepant data point does not significantly reduce 
the goodness-of-fit measure, as determined using the F-test, where xl, = /^^ a,nd 
u are the number of degrees of freedom, and = ^ [(^.(obs) — K(pred))/(Jv^(obs)] 



2 



(IBevington fc Robinsonlll992l ). 



2. The next step is to apply the requirement of contiguity: we assume that the LIC 
does not have any detached pieces that have acceptable radial velocities but cannot 
be sensibly connected to the rest of the LIC because there are lines of sight between 
the two regions that do not show radial velocities consistent with the LIC velocity 
vector. Because of the similarity of the different velocity vectors of LISM gas, coinci- 
dent projected velocities of two or more dynamical structures is common. (Those sight 
lines that have components consistent with the vector but not spatially contiguous are 
displayed as medium-sized green symbols in Figures [THIS]) . Although limiting the def- 
inition of LISM clouds to spatially coherent structures may prohibit the identification 
of complex morphologies, it has the advantage of preventing the merging of distinct 
dynamical structures with similar velocity vectors. We draw a first approximation of 
the LIC shape (see Figure [1]) consistent with all data points lying within la of the 
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predicted value if derived from high-spectral-resolution data and within 3a if drawn 
from the moderate-resolution data, with a contiguous morphology. We require at least 
4 velocity components to constitute a distinct dynamical structure. 

3. At this stage we reintroduce previously deleted velocity components that are consistent 
in velocity and continuity with the new vector. Particularly for the smallest clouds, 
some velocity components were prematurely removed from the fit at the earliest stages. 
In addition, it is at this point that other sight line properties, if available, are compared 
in order to avoid assigning clearly different collections of gas with a coincidental velocity 
vector. When several cloud vectors predict a similar projected velocity for a particular 
sight line, we assign the velocity component to the cloud with the nearest neighboring 
line of sight that is uniquely a cloud member. 

4. At this point, we have a nominal assignment of sight lines to a particular dynamical 
cloud. An iterative reevaluation of sight line membership of previously determined 
dynamical clouds is performed and occasionally a reassignment of cloud membership is 
made, although this was relatively rare. The process is then repeated for the remain- 
ing unassigned velocity components. This iterative velocity vector technique is most 
successful at identifying clouds comprised of a large number of components (e.g., LlC 
and G clouds), that subtend large angles on the sky (e.g., NGP and Mic clouds), or 
are significantly different dynamically than the average LISM flow (e.g.. Blue and Aql 
clouds). However, it has difficulty identifying compact dynamical clouds defined by 
only a handful of sight lines. In order to search for these kinds of clouds, we began the 
process with a preselected subset of sight lines which are either spatially grouped in 
a region without an identified dynamical cloud, or have a common velocity difference 



from the general LISM flow, a technique used by iFrisch et al.l (120021 ). Those that pro- 
duced a satisfactory velocity vector and survived the constraint of continuity resulted 
in the identification of some of our smallest clouds (e.g.. Dor and Oph clouds). 



2.3. Collection of Warm Nearby Interstellar Clouds 



We were able to fit rigid velocity vectors for 15 clouds in the LISMg. Absorption com- 
ponent membership and properties are given for each cloud in Tables [T]-[T5| and the velocity 
vectors and goodness -of-fit metrics, y?,, are gi ven in Table [TBI The names o f clouds are either 
historical (e.g., LIC: iMcClintock et al.lll978l : G: Ihallement fc Bertinlll992l : Blue: iGrv et al. 



^Projected and transverse velocities can be calculated for any sight line 



at 



http://cobalt.as.utexas.edu/~sredfield/LISMdynamics.htnil 
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19951 : Hyades: iRedfield fc Linskvll200ll : and North Galactic Pole fNGP): iLinskv et al.lboooh . 



or based on constellations that dominate the area of the sky coincident with the cloud lo- 
cation. The sight line members are listed in order of distance to the target star, along with 
the observed projected velocity of LISM absorption, the deviation (a) from the predicted 
projected velocity of LISM absorption, any other LISM properties along the line of sight 
(e.g., A^(H I), T, ^, D{Fe), and D{Mg)), and a list of other LISM clouds that could possibly 
explain the observed absorption component. The deviation between the observed and pre- 
dicted projected velocity is given as, a = {\vo — v^^D/cy, where Vq is the predicted projected 
velocity of absorption, is the observed velocity, and is the error in the observed veloc- 
ity, where we have imposed a minimum o"„ of 1 km s~^ for all high resolution data, and a 
minimum of 3 km s~^ for all medium resolution data. The list of other clouds that could 
possibly explain the observed absorption component were required to meet slightly relaxed 
constraints from those imposed for cloud membership, such that those clouds listed in the 
last column of Tables [THT5] are within 10° of the sight line and predict a projected velocity 
within 3(7 of the observed velocity. 

Three clouds, NGP, Oph, and Get, have > 3, indicating a relatively poor match 
between our rigid velocity vector and the observed projected velocities. We believe that the 
high xl values for these clouds indicate departures from the assumption of rigidity rather than 
the existence of several cloudlets with very similar velocity vectors. For example, the NGP 
and Oph clouds are relatively compact collections of many sight lines, 15 and 6, respectively, 
which supports a genuine connection between the absorbing material despite the poor fit to 
a rigid velocity vector. Likewise, although the Get cloud is filamentary and comprised of 
only 5 sight lines, its high velocity makes it unlikely that a set of random velocities, for a 
contiguous group of sight lines, would all be consistent with such an extreme velocity vector. 

The distributions of velocity amplitude (Vq) and direction in Galactic coordinates (/q, 
6o) for the 15 clouds are given in Figure [TBI in both the solar rest frame (heliocentric) and 
relative to the local standard of rest (LSR). The vector solutions for all clouds have similar 
directions, suggesting that there is a common history or dynamical driver for all the warm 
LISM clouds, but there i s a wide range of velo city amplitudes suggesting the presence of 
shocks in the LISM (cf. iMcGray fc Snowl Il979l ). In particular, five clouds have velocity 
components that differ significantly from the mean value: the Blue and Hyades clouds have 
Vq < 15 km s^^, and the Aql, Dor, and Get clouds have Vq > 50 km s~^. 

Figure [T7] shows the projections of the three-dimensional velocity vector solutions along 
different Galactic axes. The location of the center of each vector is placed in the direction of 
the center of the cloud at the distance of the closest star with the cloud's absorption velocity. 
This figure likewise demonstrates that the 15 velocity vectors are all variations on the same 



- 10 - 



theme, in that they all are aligned in approximately the same direction; however, significant 
differences do exist between individual velocity vectors. 

About 18.8% of the velocity components in our database cannot be assigned to any of 
our 15 derived velocity vectors. All unassigned components are listed in Table [T71 Many 
of these velocity components may represent more distant LISM clouds that subtend smaller 
fractions of the sky and are probed by too few sight lines to derive a unique velocity vector. 
Approximately 90% of the 51 unassigned absorbers are toward stars beyond 15 pc. Only five 
stars within 15 pc contain unassigned absorption components, and all have an unidentified 
absorption component in addition to absorption from identified nearby clouds. 

The first attempt t o fit a rigid velocity vector to absorption lines from nearby stars 
was made by ICrutcherl (119821 ). Seven stars that were presumed to be within 100 pc, and 
which were dominate d by a s i ngle v elocity component in moderate resolution [R ~80,000) 
Ti II observations by IStoked (Il978l ). were used to solve for a single LISM velocity vector, 
Vo = 28 km s~^, Iq = 205°, and bo = —10°. Two of the seven stars turn out to have 
Hipparcos distances >100 pc, while 4 of the remaining 5 have high- resolution UV or Ca II 
observations and are included in our database. However, all of these stars show clear evidence 



for multiple components in high-resolution spectra. Although ICrutcherl (jl982l ) was able to 
derive the general LISM flow direction, this work demonstrates that analysis of the dynamical 
structure of the LISM requires: (1) high-spectral- resolution observations of ions sensitive to 
LISM material, (2) accurate distances to the background stars, and (3) a much larger number 
of sight lines to disentangle the complicated spatial and kinematic structure of local material. 



In a series of papers, including iLallement fc BertinI (119921 ) and iLallement et al.l (119951 ). 
Lallement et al. used high-resolution {R ~1 10,000) Ca II observations and UV observations of 
Mg II and Fe II of ~16 stars to derive two rigid velocity vectors that encompassed significant 
areas of the sky. The velocity vectors were associated with the LIC, where they derive a 
solution, Vo = 25.7 km s~^, Iq = 186.1°, and bo = —16.4°, and the G Cloud, with a solution 
of Vo = 29.4 km s-\ k = 184.5°, and bo = -20.5°. The LIC and G Cloud velocity vectors 
that we derive from 5 times the number of lines of sight are similar (see Table [T6l) . For the 
LIC, we calculated a heliocentric velocity vector with a magnitude of 23.84 ib 0.90 km s~^ 
flowing toward Galactic coordinates / = 187?0 =b 3?4 and b = —13? 5 ± 3?3. Our LIC vector is 
within ~lcr of the direction proposed by ILallement fc BertinI (119921 ) but 1.8 km s~^ smaller in 



amplitude. Our G vector is almost identical to their previous determination. This agreement, 
derived from a much larger sample, demonstrates not only the reliability of the analysis, but 
also the ability to derive accurate velocity vectors from a relatively small number of sight 
lines. 



Frisch et al, 



(I2OO2I ) derived a bulk flow vector, essentially the average velocity vector 
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consistent with 96 velocity components from 60 stars within 132 pc. As we found with our 
larger database, the fit to all of the velocity components leads to a solution that approxi- 
mates the Lie velocity vector, since LIC absorption dominates observations of nea rby stars, 
although G Cloud absorption also contributes significantly. iFrisch et al.l ((2002) assumed 
the direction of that bulk flow vector for all LISM clouds (except the LIC), then identified 
compact collections of absorbers that show common velocity magnitude departures from the 
bulk flow velocity. Our Blu e Cloud vector, de rived from 10 velocity components, matches 
well with that calculated by lFrisch et al. I J2002h fro m onl y 2 velocity components. However, 
the other cloud vectors discussed by iFrisch et al.l ((2002) are not obviously comparable to 
the vectors that we derive. In particular, the remaining clouds only differ from the bulk 
flow by <3cr. Therefore, several of the velocity components identified in iFrisch et al.l (|2002l ) 
are included as members of the LIC in our calculation. Although the directions of velocity 
components in the LISM are similar, the assumption that they are identical can hinder the 
identification of distinct dynamical structures. 



In situ measurem ents derived from neutral helium (e.g^ 



Witte 



2004r). pick-up ions (e.g., 



Gloeckler et al.ll2004j ). and backscattered UV emission (e.g.. lVallerga et al.ll2004] : iLallement et al 



20041 ) , provide information on the interstellar flow vector that our solar system is presently 



encountering. iMobius et al.l (|2004l ) summarize the results of these experiments and provide 
the weighted mean values for the flow vector, Vq = 26.24ib0.45 km s~^, Iq = 183.4°±0.4°, and 
bo = —15.9° ± 0.4°. Although technically this vector is consistent within 3a of both the LIC 
and G Cloud vectors, the in situ velocity is intermediate between t he velocities of the L 
G clouds by 2.4cr and 2.9cr, respectively. Previous studies (e.g. 



Lallement fc Bertin 



C and 
1992h 



have c oncluded that the flow in the heliosphere is at the LIC velocity. iRedfield fc Linsky 
(12000! ) and others have argued that the solar system is located near the edge of the LIC 
and is moving toward the G Cloud. However, the new lower LIC velocity amplitude that we 
now derive suggests that the in situ measurements could be sampling an interaction region 
between the faster-moving G Cloud material and the slower-moving LIC cloud material, 
see Table [16] and § 15.51 This concl usion is supported b y the 6303 ± 390 K temperature 
of interstellar gas in the heliosphere (jMobius et al.ll2004j ). which, like the velocity, is inter- 
mediate between the temperature of the LIC cloud gas (7500 ± 1300 K) and the G Cloud 
gas (5500 ± 400 K), although the differences have a lower significance of 0.9cr and 1.4o", 
respectively. Additional temperature measurements, and therefore a refinement of the mean 
temperature, of both the LIC and G clouds would be possible with observations of multiple 
ions along additional sight lines. Currently, such measurer nents are available along only 29 
sight lines, 24 of which probe the LIC or G cloud material (IRedfield fc Linskyll2004bl ) . If the 
heliosphere is now located in a transition zone between the clouds, then we predict that long 
term in situ measurements will gradually approach the G Cloud velocity and temperature. 
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3. MORPHOLOGY OF WARM CLOUDS IN THE LISM 

FiguresfTHTSlshow the morphologies of each cloud encompassing the sight lines consistent 
with the 15 rigid velocity vectors derived from the LISM databasqf|. The sight lines utilized 
in the velocity vector calculations are indicated by the large blue symbols, while sight lines 
with consistent projected velocities, but assigned to other clouds, are shown by the medium 
green symbols, and those sight lines with observed projected velocities that are inconsistent 
with the velocity vector fit are indicated by the small red symbols. The projected morphology 
of each cloud is drawn to include all sight lines used in the velocity fit, while avoiding all 
lines of sight that are inconsistent with the velocity vector. Each figure shows the cloud 
morphology from four different directions in Galactic coordinates. 

Although the filamentary nature of the clouds could be exaggerated in some cases be- 
cause of the low spatial sampling and our requirement of cohesion, approximately a third 
of the clouds have projected morphologies that are clearly filamentary. It is possible that 
a couple of the "compact" morphologies may actually be filamentary, but due to a chance 
orientation, are projected as a compact cloud on the sky. However, it would be highly un- 
likely that the orientations of many "compact" clouds would be precisely aligned along the 
line of sight to hide the true morphology of the clouds. Determining the true morphology 
of these clouds, regardless of orientation, requires a database with high spatial and distance 
sampling. The orientations of the observed filamentary clouds are not similar, which argues 
against an association with a global magnetic field that may thread through the LISM. In- 
stead, the filamentary regions, which generally trace the boundary between the LIC and G 
clouds, may indicate regions of cloud-cloud interactions, where the rigid velocity structure 
is disrupted and potentially shocked by the collision of two adjacent clouds. 



4. PHYSICAL PROPERTIES 

Table [T8] summarizes basic physical properties (e.g., coordinates of the cloud center, up- 
per limits to the cloud's distance, projected surface area, weighted mean temperature, turbu- 
lent velocity, and depletion of iron and magnesium) of the 15 clouds. The surface area on the 
sky is simply the surface area of the projected boundaries shown in Figures [THTH The temper- 
ature and turbulent velocities are derived from comparisons of the measured Doppler widths 
of absorption lines of elements with different atomic masses (e.g., deuterium and iron) using 



■^Probable cloud membership based on the projected cloud morphologies, can be calculated for any sight 
line at http://cobalt.as.utexas.edu/~sredfield/LISMdynamics.html , 
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the relation, 6^ = 2kT/m + ^'^ flRedfield fc Linskvll2004bl : IWood et al.lll996h. The depletions, 



D(JQ = ^og{X/H) — log (X/if)^, are calculated using the lAsplund. Grevesse. &: Sauval 
(I2OO5I ) solar abundances, where the hydrogen abundance is typically calculated from D I 
and converted to H I using the remarkably constant LISM D/H ratio of 1.56 ± 0.04 x 10~^ 
( iLinsky et al.ll2006l ). The depletions do not take into account partial ionization of hydrogen, 
which is likely important, or neutral or doubly ionized magnesium and iron, which are likely 
much less importar it since they are not expected to become a dominant ionization stage of 
either element (cf. ISlavin fc Frischl l2002l : iLehner et al.ll2003l ). The weighted mean and la 
uncertainty of the mean are listed for all physical properties with measurements on multiple 
sight lines. No significant correlation appears to exist between any of the physical properties 
listed and cloud morphology. In the following discussion of physical properties of individual 
clouds, we consider only the nine clouds that have more than 1 sight line with physical 
measurements. 



4.1. Distance Limits 

Although a detailed look at the distances of the 15 LISM clouds is beyond the scope 
of this paper, we can immediately place distance constraints based on the distance of our 
background sources and provide some insight into the three dimensional structure of the 
LISM. In Table [TSl we list the distance of the closest star with absorption from the cloud, 
which provides an upper limit to the distance of the cloud. All of the clouds lie within 15 pc, 
and half lie within ~5 pc, which is much smaller than the volume of the Local Bubble, but 
consistent with the large projected surface area that these clouds subtend. The distribution 
of many clouds with a range of dynamical properties in such a small volume, makes collisions 
between clouds a real possibility. The implications of such interactions are explored in § O 

Although the G Cloud has a more stringent distance upper limit than the LIC, since 
the temperature and velocity of the interstellar mate rial that is flowing i nto the solar system 



is consistent with early estimates of LIC material (iMobius et al.ll2004l ). implying that the 
Sun is currently inside the LIC (although, see § 15.51 for a detailed discussion of this topic). 
However, since LIC absorption is not observed in all directions (e.g., toward the Galactic 
Center) and since the Sun is moving in roughly the direction of the Galactic Center and the 
G Cloud (Figure [T^. th e heliosphere has been thought to be at the very edge of the LIC 



(IRedfleld fc Linskvll2000l ). 
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4.2. Volume Filling Factor of the Warm LISM 



We have not yet created a full three-dimensional morphological model of the LISM, 
but with a few assumptions and a simple toy model, we can estimate the volume filling 
factor of the warm partially ionized gas in the LISM. First, we assume that all of the warm 
LISM material is located within 15 pc of the Sun. Although the LISM is often considered 
to be the volume of material within the Local Bubble, which extends out to roughly 100 pc 
in all directions, it seems that most of the warm materi al is located only a short d istance 
from the Sun. This is shown, for example, in Figure 14 in lRedfield fc LinskyI (j2004al ). where 
the average number of absorbers per unit distance levels off at a distance of ~15 pc. No 
significant correlation exists between the observed line width or column density and the 
distance of the background star. Therefore, it is unlikely that unrecognized line blends along 
more distant sight lines are the cause of the observed leveling off of the average number 
of absorbers, but is indicative of the t rue distribution of warm gas in the LISM. Based on 
Figure 14 of iRedfield &: LinskyI (j2004a|), there are on average ~1.7 absorbers per sight line. 
Therefore, the projected surface area of all LISM clouds should total ~1.7x47r. 

Initially, we assume that all of the warm LISM clouds are similar in size and and density 
to the Lie. We assume that all warm clouds have a radius of 1.5 pc and a total hydrogen 
density of 0.2 cm~^, obtai ned from measuremen ts of the He I volume density streaming 
into the solar system from iGloeckler et al.l (120041 ). the H I to He I ratio of the LISM from 
Dupuis et al.l (Il995l ). and the three-dimensional model of the LIC from iRedfield &: Linsky 
( 120001 ). In addition, the assumed radius and density result in a full cloud hydrogen column 
density of ~2xl0^^ cm~^ , which matches the typical observed hydrogen column density 
jRedfield fc Linskvll2004j ). 



We ran 1000 simulations of cloud distributions, for a range of total warm clouds (n) 
in the LISM from n = 1-100 clouds, where all clouds were randomly centered at distances 
from 0-15 pc. Since clouds are not allowed to overlap, only one cloud can surround the 
solar system, that is, a LIC analog. In order to more closely match the model with the 
observed LISM, we have assumed the LIC projected surface area of 18270 square degrees 
for the LIC analog, see Table [18], instead of 41250 square degrees. The remaining clouds 
projected surface areas were calculated based on their size, geometry, and distance. In the 
spherical cloud scenario, the solution to an average of 1.7 absorbers per sight line leads to 
~55 clouds within 15 pc, and a volume filling factor of warm partially ionized material of 
-5.5%. 

We also ran a suite of simulations varying the geometry (i.e., ellipsoids with a range of 
aspect ratios from 1.33:1 to 10:1, both flattened (i.e., pancakes) and elongated (i.e., cigars)) 
and fraction of ellipsoid to spherical (which ranged from 0.3 to 1.0). The orientation of all 
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clouds were determined randomly. Non-spherical geometries naturally lead to larger volume 
filling factors and fewer clouds, since a greater total volume can be produced with fewer 
clouds without necessarily increasing the projected surface area. The volume filling factors 
that resulted ranged from 5.5% to 19%. 

Figure [T8] compares the projected surface area distribution of our sample with two of our 
idealized simulations. The distribution of observed projected surface areas matches fairly 
well with that predicted for both the ~55 spherical LIC-like clouds within 15 pc, and a 
simulation of 35 clouds, half of which are ellipsoids with an elongated aspect ratio of 10:1. 
There is an observational bias toward detecting the nearest clouds with the largest projected 
surface areas. We restricted our dynamical cloud modeling to collections of gas that had at 
least 4 sight lines with which to determine a velocity vector. 

About 18.8% of absorbers are not accounted for in our 15 cloud dynamical model of 
the LISM. The missing absorbers may represent detections of more distant and smaller 
projected surface area clouds. If we assume that we have detected all clouds with a log 
surface area > 3.1 (~1260 square degrees), which is the lower limit of the LISM clouds with 
measured velocity vectors (see Table [T5]) . we can estimate the number of sight lines that 
probe "undetected" clouds, or clouds with a log surface area < 3.1, (see Figure [T8l). Our toy 
models have 10.2-11.9 clouds with log surface area > 3.1, slightly lower, but similar to the 
15 observed clouds. The total projected surface area for all clouds with log surface area > 
3.1 range from 55290 to 60430 square degrees in our simulations. This matches well the total 
projected surface area of the observed clouds (57830 square degrees). With an estimate of 
the total projected surface area of "undetected" clouds, and the assumption that all sight 
lines are uniformly distributed, we can predict the percentage of observed components that 
will be left over, after those associated with the large nearby clouds are removed. The 
percentage of unassociated velocity components in the simulations are between 14.5% and 
21.5%, which matches closely the percentage of components in our database (18.8%) that 
are unassigned, which allows for the possibility that the absorbing material associated with 
these components is indeed located within ~15pc, even though the background star is much 
further away. 

5. WHERE CLOUDS COLLIDE 

5.1. The "Ring of Fire" Around the G Cloud 

Figure [11] shows the projected morphologies of all 15 clouds, which collectively cover 
more than 90% of the sky. The LIC and G Clouds clearly dominate the sky, and contain 
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large areas where one or the other is the only absorber along the line of sight. At the 
boundaries of the LIC and G clouds, several overlapping absorbers are typically present, 
particularly near / from 40° to 80° and b from -15° to +30°, as well as near / from 270° 
to 320° and b from +20° to +50° and from -70° to -30°. These areas at the boundaries 
of the LIC and G clouds may be dynamical interaction zones, which produce "new" clouds 
with significantly different kinematic properties. We refer to the active boundary of the G 
Cloud where the G and LIC clouds may be colliding as the "Ring of Fire" , in analogy to 
the Pacific Ocean "Ring of Fire" where dominant tectonic plates (here interstellar clouds) 
interact, resulting in a highly dynamic interaction zone that gives rise to earthquakes and 
volcanos (here interstellar shocks, heating, or turbulent flows). 

One example of this interaction may be the Mic Cloud whose morphology appears to 
mirror the adjacent sections of the LIC and G clouds, as shown in Figure [201 At positive 
Galactic latitudes, where the projected morphologies of the LIC, G, and Mic clouds are 
coincident, the median predicted radial velocity difference between the LIC and G clouds is 
~5.5 km s~^. The Mic Cloud may have been created by the faster G Cloud colliding with 
the LIC, which is moving ~5.5 km s~^ slower in the radial direction. 



5.2. Cloud Interactions, Turbulence, and Shocks 



Except for the LIC, many clouds have only 3-5 sight lines with measurements of physical 
properties, and six of the clouds have one or no sight lines with measured physical properties. 
As a result, it is difficult to explore the homogeneity or variation of properties across an 
individual cloud. Small-scale variations are not observed in the warm LISM clouds, based on 
identical absorption proper ties of nearby binary stars (e.g., a Gen A and B.lLinskv fc Wood 



1996: Lallement et al. 1995; and a CMa A and B Lallement et al. 1994: Hebrard et al 



19991 



and by the lack of significant var iation among a sample o f 18 closely spaced Hyades stars, 
down to scales between 0.1-1 pc (iRedfield &: Linskyll200ll ). Therefore, we may expect that 
the physical properties of LISM clouds are relatively homogeneous, or at least slowly varying 
within a cloud. 

Table [18] lists the weighted mean values of physical properties for all cloud members, as 
well as the weighted average standard deviation, which gives an indication of how tightly the 
values are scattered about the weighted mean. For example, D(Fe) in the LIC, even though 
there are 12 measurements, much more than any other cloud, it has the lowest weighted aver- 
age standard deviation, indicating that little variation is detected across the LIC. In contrast, 
due to the a wide range of values in the Hyades and Mic clouds, there is a large weighted 
average variance (although these clouds have only a few measurements). In particular, note 



-17- 



the anomalous depletion measurement of G191-B2B associated with the Hyades Cloud, see 
Table [71 The ~8.6 km s~^ component, consistent dynamically and spatially w ith the Hyades 



Cloud, is detected in both low- ionization ions (e.g., D I, N I, O I, Mg II. etc. iLemoine et al. 



19961 : iRedfield fc LinskyI l2004al ) and high-ionization ions (e.g., C IV, IVennes &: Lang|200ll ). 
The nature of the absorbing material along this line of sight is not well known, and the 
high-ionization material may be asso ciated with nebular circumstellar material surround- 



ing G191-B2B (IBannister et al.l 120031 ). Although such contamination may be present along 
some sight lines in our sample, the need in this analysis to bring together a large number of 
independent LISM measurements aids in reducing and identifying anomalous data points. 
Indeed, the high weighted average standard deviation of depletion in the Hyades Cloud, 
clearly identifies the G191-B2B sight line as anomalous. Henceforth, we assume all measure- 
ments sample the physical properties of the LISM, although highly deviant data points may 
indicate interesting sight lines that require additional observations and further attention. 

If we assume we can estimate an individual cloud's mean properties by assuming that 
clouds are approximately homogeneous, then we can use the few available physical mea- 
surements to make a reasonable estimate of the cloud properties, and search for possible 
correlations. For example, the Mic Cloud, already identified by its filamentary morphol- 
ogy and location at the boundary of the LIC and G clouds, has the highest temperature 
((T) = 9900 K) and one of the highest turbulent velocities ((^) = 3.1 km s~^). These prop- 
erties support the argument that the Mic Cloud is the result of the collision of the LIC and 
G Clouds. 

The weighted mean depletions of iron [Z}(Fe)] and magnesium [D(Mg)] are plotted versus 
the weighted mean turbulent velocity (^) for the nine clouds in Figure [211 A clear correlation 
is evident. For iron, the linear correlation coefficient, r = 0.69, and the probability (Pc) that 
this distribution could be drawn from an uncorrelated parent population is 1.7%, while for 
magnesium, r = 0.73 and Pc = 1.2%. 

The correlation of small depletions with high turbulence suggests that the destruction 
of dust grains has returned these ions to the gas phase. A possible alternative explanation is 
that statistically, clouds with higher turbulence have higher percentage ionization of hydro- 
gen, since the depletions were computed assu ming that hydrogen is n eutral. However, using 



the ionization model of the LIC produced by [Slavin fc FrischI (|2002[ ). taking the ionization 



of hydrogen, magnesium, and iron into account only produces a 0.05 to 0.10 decrease in the 
measured depletion. This adjustment is significantly less than the typical la error for the 
weighted mean depletion for individual clouds, and much less than the ~1 dex variation 
seen over all LISM clouds. Therefore, since we currently have no evidence for a correla- 
tion between turbulence and ionization structure and the depletion adjustment using LIC 
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ionization models is small, we consider here possible dust destruction explanations for the 
observed correlation. Dust destruction in the w arm partially ionized ISM i s often discussed i n 



the context of shocks produced by supernovae (ISavage fc Sembachlll996l : I Jones et al.lll994j ). 
Shocks may also be produced in the LISM from either turbulent motions possibly driven 
by shear flow interactions between clouds, or from direct macroscopic collision s of cl ouds. 



The thermal sound speed Cs = \ /nkT/p, using LISM dens ities from iRedfieldl ( 20061 ) and 



the mean LISM temperature from iRedfield fc LinskyI (j2004bl ) , is ~8 km s ^ . For the clouds 



with the highest turbulence (e.g., Mic and Eri), the sight line-averaged turbulent velocity is 
~3.5 km s^^, which results in a turbulent Mach number (Mg) of ~0.4. Although there may 
be regions of enhanced turbulent motions, perhaps at the interaction boundaries of clouds, 
average turbulent velocities are not high enough to produce shocks. 

Macroscopic velocity differences between the 15 LISM clouds can be significantly greater 
than the sound speed if they are interacting. Figure [22] shows the distribution of predicted 
velocity differences (AV^) between the 15 LISM clouds when multiple clouds are predicted to 
lie along a line of sight. For a uniform sample of hypothetical individual sight lines across the 
entire sky, we predict the number of the 15 LISM clouds are predicted to lie along the line 
of sight based on their boundaries shown in Figures [THTBl If multiple clouds are predicted 
to lie along the line of sight, we calculate all possible cloud velocity differences, which are 
shown in Figure [221 For example, a hypothetical sight line that traverses two clouds (e.g., 
at / = 270° and b = 0° where the G and Get clouds overlap) will provide one velocity 
difference measurement in Figure [251 whereas, if three clouds are predicted along the line of 
sight (e.g., at I = 165° and 6 = 0° where the LIG, Aur and Hyades clouds overlap), three 
possible velocity difference measurements are shown. Galculating the radial and transverse 
components of a velocity vector of LISM material along an arbitrary line of sight requires 
projecting the velocity vector (Vq, Iq, bo) along the radial and transverse unit vectors in the 
arbitrary direction (/, b). The magnitude of the radial and transverse velocities are calculated 
from 

Vr = Vo(cos6cos6ocos(/o — + sin 6o sin 6), (1) 
Vi = Vocosbosm{lo - I), (2) 
Vb = Vo(sin6ocos6 — cos 6o sin 6 cos(/o — 0)- (3) 

Since we do not have a fully three-dimensional model of these LISM clouds (i.e., we do not 
know which clouds are in fact adjacent and actively interacting), we do not know whether all 
these velocity differences are actually realized. However, a large percentage of the possible 
velocity differences are supersonic and thus could cause shocks where adjacent clouds meet. 

The question of grain destruction/erosion involves many factors including grain size, 
composition, porosity, the relative speed of collisions with other grains or particles, com- 
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pression ratio in shocks, grain charge, ma gnetic fields, and turbulence. This to pic has been 
addressed by a number of authors (e.g., I Jones. Tielens. fc HoUenbachl Il996l ). iFrisch et al. 



( 119991 ) and others showed that the observed Fe and Mg ions in warm gas like the LIC comes 
from grain destruction whether by shocks, grain-grain collisions, or other physical processes. 
Since most of the Fe and Mg in the LIC is locked up in grains, small differences in grain 
destruction between or within clouds can produce large differences in the gas phase abun- 
dances of these elements. High-speed supernova-generated shocks (50-200 km s~^) are often 
cited as the main grain destruction method, but grain-grain collisions with relative velocities 
exceeding only 2.7 km s~^ for silicate grains or 1.2 km s~^ for carbonaceous grains can lead 
to grain shattering (jjones et al.lll996l ). Velocities exceeding these values are typically found 
between two clouds along the same line of sight (Figure [22]) and are similar to the measured 
cloud turbulent velocities. Given t hat interstellar dust grain s are t ypically charged and the 
ISM is magnetized and turbulent, lYan. Lazarian. &: Draind (120041 ) showed that MHD tur- 
bulence can accelerate the grains through gyroresonance interactions leading to supersonic 
grain speeds, grain-grain collisions and shattering. This process could be the physical basis 
for the observed correlation of low metal depletions with high turbulent velocities. 



5.3. Connection with Radio Scintillation Screens 



For many years, radio observers have called attention to a rapid variability of certain 
quasars and pulsars on hourly-to-yearly timescales that has been attributed to interstellar 
scintillation. The scattering screens responsible for the scintillation are generally assumed 
to be turbulent regions of enhanced electron density. Extensive monitoring of a source and 
measurement of time delays as seen by widely separated radio telescopes provide critical data 
for estimating the distance to the scattering screen, as well as its size, transverse velocity, and 



shape . In their VLA survey of northern sky AGN for rapid intraday variability, iLovell et al. 



( 120071 ) found that 56% of the sources are variable on timescales of hours to several days, 
but rapid variability indicative of nearby scattering screens is rare indicating that nearby 
scattering screens cover only a small fraction of the sky. 

Studies of intraday variability of three quasars (J1819+3845, PKS 1257-326, and PKS 
0405-385) and two pulsars (PSR J0437-4715 and PSR B1133-M6) find that some scattering 
screens lie within the LISM at short distances from the Sun, although the distances have 
signifi cant uncertainty and are model de pendent. For example, iDennett-Thorpe fc de Bruyn 
(2003) and Macquart &: de Bruvn (I2OO6) esti r aate t hat the scattering screen toward J1819-I-3845 
is only 1-12 pc from the Sun. iBignall et~aL (2006) found t hat the scattering screen toward 



PKS 1257-3826 hes at a distance somewhat closer than 10 pc. iRickett. Kedziora-Chudczer. fc Jauncey 
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( 120021 ) place the anisotropic scattering screen toward PKS 0405-385 at between 2 and 30 pc 
from the Sun with a pref erred distance of 25 pc, though an alysis of more recent data by 
Kedziora-Chudczerl (120061 ) suggests a distance of about 10 pc. ISmirnova. Gwinn. fc Shishov 
( 120061 ) show that the scattering screen toward PSR J0437-4715 also lies at about 10 pc from 
the Sun and is likely the same screen that causes the scintillation of PKS 0405-385, which 
is only 10 degrees away. If so , this scattering screen is extended rather than very restricted 



m size. iPutney fc Stinebrind (120061 ) present six pulsars that show evidence of multiple scin- 
tillation screens along their line of sight, the vast majority of which are located well beyond 
the Local Bubble. One of their nearest pulsars, PSR B1133+16, shows evidence for a nearby 
scintillation screen only 21.6 pc from the Sun. Several other scintillating quasars, or intraday 
variables, show annual cycles tha t may provide cons t raints on the distance to th e scintilla- 
tion scree n, such as B0917-F6 24 jRickett et allboOll : Ijauncev fc MacquartI [200ll ) and PKS 
1519-273 Jjauncev et al.ll2003h . 



We find that the five nearby scattering screens all lie close to the edges of several of 
our dynamical clouds, as indicated in Figure [19], where the direction of the radio scintil- 
lation sources are indicated by star symbols. In particular, three of the five lie near the 
interface of the LIC and G clouds. The radial velocity differences between the LIC and G 
clouds in these directions are generally quite small (i.e., ~1 km s~^). However, the trans- 
verse motion differences between the LIC and G clouds in these directions can be quite 
substantial, reaching 6-7 km s~^. These regions of significant transverse velocity differences 
could induce shear flows and generate turbulence. The annual variation of the scintillation 
timescale of intraday variables, is a function of the diffraction pattern of the screen and its 
transverse velocity. With our rigid velocity vectors of LISM clouds, we are able to calculate 
the transverse motions of local clouds. Even with just a handful of sight lines, we can in- 
vestigate the relationship between the LISM and the scintillation screens, but many more 
radio scintillation and high resolution LISM absorption line observations are needed together 
with a fully three dimensional morphological and kinematic model of the warm LISM, in 
order to fully explore the ph y sical connection between warm clouds and scintillation screens. 
Linsky. Rickett. fc Redfieldl (120071 ) more fully explore the relationship between scattering 
screens and LISM clouds, including a direct comparison of the transverse velocities of the 
screens and the clouds. 



5.4. Connection with Cold Dense Structures in the LISM 



In their 21-cni abs orption line study of the warm and cold neutral interstellar gas, 
Heiles &: Trolandl (120031 ) mapped a region of cold gas centered at I = 225°, b = 44° that 
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extends over 30° in Galactic longitude. They found that the gas temperature is about 25 K. 



Meyer et al.l (120061 ) observed narrow Na I absorption due to this cold cloud toward a series of 
nearby stars and confirm a cloud temperature ~20 K and turbulent velocity of ~0.4 km s~^. 
Based on the well-known distances of the observed stars, they were able to show that the 
distance to this cold gas must be less than 41 pc, with a corresponding aspect ratio (length 
perpendicular to the line of sight versus length along the line of sight) of 70:1. However, the 
cloud could be as close as 2 pc. Thus the cold gas structure is located inside the LISM. 

The Galactic coordinates of the cold gas correspond to a region that is not clearly a 
part of any individual warm cloud (see Figure [T9|) . but near the boundaries of several clouds, 
including the LIC, G, Aur, Gem, and Leo clouds. In particular, the Gem Cloud has a high 
radial velocity in the direction of the cold cloud, ~24 km s~^, whereas the other clouds 
have modest radial velocities, ranging from 8-12 km s~^. The resulting high radial velocity 
differences are indicated in Figure [221 Along this line of sight, the high radial velocity (~24 
km s~^) Gem cloud may be compressing material as it collides with slow moving Leo and 
Aur (~12 km s~^), and ultimately the L IC (~10 km s ~^). T he heliocentric velocity (~11.5 



km s ^) of the cold material observed by lMeyer et al.l (l2006l ) matches well with the velocity 



of the slow moving Leo and Aur clouds, as expected if it was formed by the compression 
of the Gem cloud against the Leo and Aur clouds, and the cold material may actually be 
physically associated with the warm material observed in the Leo and Aur clouds. 

The collision of warm gas clouds to produce small sheetlike cold neutral clouds has been 



explored through detailed simula tions of a turbulent interste ll ar me dium (e.g.. lVazquez-Semadeni et al 



20061 : lAudit fc Hennebel"l3 boosh . Also iMcKee fc Ostrikeil Jl977h predicted that cold 



neu- 



tral clouds must be surrounded by warm clouds (in pressure equilibrium in their model) to 
shield the cold gas from UV and X-ray heating and ionization. Our rigid velocity vector 
solutions certainly indicate that in the direction of the cold cloud, relatively large kinematic 
differences exist between clouds (e.g, the Gem cloud is moving ~12 km s~^ faster than the 
Leo and Aur clouds, and ~14 km s~^ faster than the LIC and G clouds). Given that the 
farthest distance limit for the five clouds in the proximity of the cold cloud line of sight is 
11.1 pc (the Leo Cloud), and that the limit on the Gem Cloud is 6.7 pc, these five clouds 
are in very close proximity in distance as well, and collisions between these warm clouds are 
likely. In particular, it is critical to have a collision of material along the radial direction to 
maximize the chances of detection. Because of the extreme aspect ratio of this cloud, if it 
were oriented along the line of sight, the projection on the sky would be extremely small. 
Further work on the distances of both the warm LISM clouds and the cold cloud is needed 
to determine whether there is a spatial and dynamical connection between these interstellar 
structures. 
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5.5. The Transition between the LIC and G Clouds 

Until now, we have assumed that individual clouds are rigid structures, each with a 
simple velocity vector characterizing all of the included sight lines within the radial velocity 
measurement errors. This simple approximation may not be valid, as shown by two sets of 
data. One, noted in Section 12. 3^ is that the velocity and, to a lesser degree of significance, 
the temperature of interstellar gas flowing through the heliosphere are intermediate in value 
between the corresponding quantities in the LIC and G clouds (see Table [T6|) . implying that 
the heliosphere lies in a transition zon e between the two c louds where there is a gradient in 



properties. The other evidence is that iRedfield fc LinskyI (l200ll ) noted that the components 



assigned to the LIC ir i the direction of the Hyad es have radial velocities 2.9 km s ^ smaller 



than predicted by the iLallement fc BertinI (119921 ) LIC vector. These smaller radial velocities 



suggest a deceleration of the LIC flow in the forward direction, where it may be interacting 
with the Hyades Cloud. We now find that these absorption components have radial velocities 
~1.0 km s~^ smaller than predicted by the new LIC vector. As a test we removed these 16 
components from the LIC vector calculation and found that the vector velocity amplitude 
increased by only 0.39 km s^^. This does not change our conclusion that the heliosphere 
is in a transition zone between the LIC and G clouds, or that the LIC is decelerated at its 
forward edge. 

We also considered whether the LIC and G Clouds are really one cloud with a gradient 
of physical properties across their combined length. We tested this hypothesis by plotting 
the physical parameters for the LIC and G cloud sight lines with respect to angle relative 
to the downwind direction (Figure [231) and with respect to the hydrogen column density 
(Figure m]). No correlation exists between angle and hydrogen column density. Since the 
LIC Cloud is mostly in the downwind direction and the G Cloud mostly in the upwind 
direction, the angle from the downwind direction is a discriminant between the two clouds. 
With only one exception, the gas temperatures for LIC sight lines are all larger than for the 
G Cloud sight lines, implying that the two cloud approximation is valid. We note, however, 
a trend of higher LIC Cloud temperatures with increasing A^(II I) and toward the cross- 
wind direction where the two Clouds meet. The turbulent velocities do not show a definite 
trend with angle or A^(II I). If there were a velocity gradient through these clouds, one 
would expect larger line broadening (and thus higher turbulent velocity) in the downwind 
and upwind directions, which is not seen. Finally, the metal depletions, [Z)(Mg) and D{Fe)], 
show clear trends of decreasing with larger A^(H I) and increasing with angle, with both the 
LIC and G Cloud sight lines fitting these trends. On average, the LIC is significantly more 
depleted than the G cloud. The trend of decreasing gas-phase abu ndances as a function of 



increasing N(B. I) is well docume nted along distant sight lines (e.g.. IWakker fc Mathisll2000 



Jenkins. Savage, fc Spitzerlll986l ). This implies that the processes that remove and replace 
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ions to the gas phase occur on scales smaller than the LIC and G clouds, and that the volume 
density of the LIC and G clouds may not be constant. Alternatively, decreasing D{Mg) and 
D(Fe) with increasing iV(H I) may be explained if H I is photoionized at the edges of clouds 
and is increasingly neutral with increasing H I through the centers of clouds due to shielding 
of UV radiation from the cloud itself. Thus, while the fraction of magnesium and iron in the 
gas phase relative to the total amount of hydrogen remains the same throughout, D{Mg) 
and D{Fe) will decrease as hydrogen becomes predominately neutral. On the basis of the 
temperature and turbulent velocity data, we conclude that the evidence supports the idea 
that the LIC and G Clouds are separate entities with their own distinct properties, but 
there is likely a narrow transition zone between the two clouds where the heliosphere is now 
located. 



6. CONCLUSIONS 



We have created a database consisting of interstellar radial velocities and gas physical 
properties for 157 sight lines toward stars within 100 pc of the Sun. The data were extracted 
from high-resolution UV spectra obtained with the GHRS and STIS instruments on the HST 
and ground-based Ca II spectra. This database has allowed us to create a dynamical model 
of the local interstellar medium0 including 15 warm gas clouds, which we define as contiguous 
parcels of interstellar gas with consistent kinematical properties. Although measurements 
of physical properties are sparse, for the LIC, which has the most such measurements, the 
properties seem to be homogeneous. Using this database, we find that: 



1. The flow velocity vectors for these 15 clouds fit 81.2% of the velocity components in 
the database to within the radial velocity measurement errors. These clouds all lie 
within 15 pc of the Sun. The remaining velocity components may be produced in more 
distant clouds that subtend smaller angles with less than the four lines of sight needed 
to compute a useful velocity vector. 

2. The directions of most of these velocity vectors are roughly parallel with their flow 
from the Scorpio-Centaurus Association. The velocity amplitudes have a considerable 
range, leading us to compute relative velocities between adjacent clouds that are often 
supersonic and therefore capable of producing shocks. 



^Projected and transverse velocities, along with probable cloud membership based 
on the projected cloud morphologies, can be calculated for any sight line at 
http: / / cobalt.as.utexas.edu/~sredfield/LISMdynamics.html 
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3. About one-third of the clouds appear to have filamentary structures. 

4. All of the clouds for which we have physical properties along three or more sight 
lines are warm with mean temperatures in the range of 5300-9900 K, although the 
uncertainties in these measurements are often large. We estimate that between 5.5% 
and 19% of the LISM within 15 pc of the Sun is filled with warm gas clouds. 

5. We find a strong correlation of low metal depletion with large turbulent velocity. Since 
high turbulence suggests the presence or recent existence of shocks, this correlation 
could be explained by shock dissipation of dust grains that returns the metals to the 
gas phase. 

6. Contrary to previous work, the heliosphere appears to be located in a transition zone 
between the LIC and G Clouds. The evidence for this is that the temperature and 
velocity of the interstellar gas flowing through the heliosphere are both intermediate 
between these quantities measured in the LIC and G Clouds. The deviation in veloc- 
ity ranges from 2.4-2.9cr and in temperature the deviation is less significant ranging 
from 0.9-1.4(7. Additional observations of multiple ions are required to increase the 
number of temperature measurements of the LIC and G clouds in order to increase the 
significance of any possible deviation between these clouds and in situ measurements. 
Previous work based on much smaller velocity data sets placed the heliosphere inside 
but near the edge of the LIC. 

7. The G Cloud is surrounded by and likely interacting with a number of other clouds. We 
refer to this active boundary as the "Ring of Fire" . The filamentary-shaped Mic Cloud 
has the same shape as the boundary of the G and LIC clouds and may be indicative 
of an interaction between these two clouds. The high temperature and turbulence of 
the Mic Cloud support this conclusion. 

8. The nearby scintillation screens toward three quasars and two pulsars are located near 
cloud boundaries, and three of the five are in directions where the LIC and G Clouds 
may be interacting. The large transverse relative velocities between these two clouds 
could produce the turbulence that is the cause of the scintillation. 



9. The ri earby cold cloud recently observed by iHeiles fc Trolandl (120031 ) and iMeyer et al. 



( 120061 ) is in a direction where it could be surrounded by several warm clouds. We find 
evidence for significant compression based on large macroscopic velocity differences 
between warm clouds in the direction of the cold cloud. The alignment of the cold 
cloud matches well with the alignment of the high-velocity Gem Cloud, which may 
be colliding with the slower moving Leo, Aur, and LIC clouds. Compression of warm 
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material may be the origin mechanism for such an isolated cold cloud in the Local 
Bubble. 
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Fig. 1. — Four projections of the LIC boundaries in Galactic coordinates. All sight lines used 
in our analysis are displayed. The large blue symbols indicate sight lines used to calculate 
the LIC velocity vector. The medium- sized green symbols indicate sight lines with projected 
velocities that are consistent with the estimated vector, but are considered part of another 
cloud (see § 2.2), while the small red symbols indicate lines of sight that are inconsistent 
with the calculated velocity vector. The boundaries of the LIC are drawn to encompass all 
consistent sight lines (i.e., blue symbols), while avoiding all other lines of sight (i.e., red and 
green symbols) . The upwind heliocentric direction of the velocity vector is indicated by 
the ® symbol, while the downwind heliocentric direction is indicated by the symbol. The 
four projections from upper left and moving counter-clockwise are: a Hammer projection of 
Galactic coordinates with the Galactic center in the middle, a Hammer projection of Galactic 
coordinates with the Galactic anti-center in the middle, a Lambert projection from the south 
Galactic pole, and a Lambert projection from the north Galactic pole. 
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Fig. 2. — Same as Figured] but for the G Cloud. 
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Fig. 3. — Same as Figured] but for the Blue Cloud. 
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Fig. 4. — Same as Figure [T] but for the Aql Cloud. Note the small number of coincident 
velocities (i.e., green medium-sized symbols). The Aql Cloud velocity vector is significantly 
different than the average LISM flow, but successfully characterizes 9 closely spaced sight 
lines. 
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Fig. 5. — Same as Figure [T] but for the Eri Cloud. 
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Fig. 6. — Same as Figured] but for the Aur Cloud. 
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Fig. 7. — Same as Figure [H but for the Hyades Cloud. 



- 38 - 




Fig. 8. — Same as Figure [H but for the Mic Cloud. 
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Fig. 9. — Same as Figure [T] but for the Oph Cloud. 
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Fig. 10.— 



Same as Figure [T] but for the Gem Cloud. 
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Fig. 11. — Same as Figure [T] but for the NGP Cloud. 
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Fig. 12.— 



Same as Figure [U but for the Leo Cloud. 




Fig. 13. — Same as Figure [U but for the Dor Cloud. 




Fig. 14. — Same as Figure [T] but for the Vel Cloud. 
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Fig. 15. — Same as Figure [T] but for the Get Cloud. Note that the boundaries of the Get 
Cloud include a couple sight lines that were not used in the velocity vector calculation which 
nonetheless have consistent projected velocities (i.e., green symbol sight lines), and therefore 
may traverse Get Cloud material. 
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Downwind LISM Velocity Vectors Relative to the Local Standard of Rest 
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Fig. 16. — Distributions of fit parameters, downwind heliocentric [top) velocity (Vo), and di- 
rection in Galactic coordinates [Iq, bo) for 15 clouds identified within 15 pc. The distribution 
of downwind velocity and direction relative to the Local Standard of Rest (LSR) is shown in 
the bottom panels, where the upstream solar motion relative to t he LS R {Vq = 13.4 km s~^; 
^0 = 207.7°; 6q = —32.4°), was derived by iDehnen &: Binneyl (119981 ). The bin sizes in Vq 
are 4 km s~^, 9° in Iq, and 6° in fog- AH velocity vectors appear to be driven in the same 
direction, although at a range of velocities. 
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Fig. 17. — Heliocentric velocity vectors of all 15 clouds. The vectors are centered in the 
direction of the center of the cloud and at the distance of the closest star with the cloud's 
absorption velocity and point downwind. The Sun is moving in roughly the opposite direction 
as the LISM clouds. Starting from the top left and moving clockwise, the plots are viewed 
from / = 230° and 6 = 45°, the North Galactic Pole (6 = 90°), the Galactic Center (/ = 0° 
and 6 = 0°), and Galactic East (/ = 90° and 6 = 0°). 



-48- 




1x10^ 2x10^ 3x10' 

Surface Area (sq. deg) 



3.0 3.5 4.0 

log Surface Area {sq. deg) 



Fig. 18. — Histogram (black) of the observed angular areas of the 15 nearby dynamical 
clouds, based on their projected morphologies shown in Figures flHTSl The histograms on 
the left have bin sizes of 1000 square degrees, while the histograms on the right have bin sizes 
of 0.1 dex in logarithmic square degrees. The red histogram indicates the average distribution 
of angular areas of the simple model of 55 randomly distributed spherical LIC-like clouds 
within 15 pc of the Sun, discussed in Section This model leads to a volume filling factor 
of ~5.5%. The blue histogram shows another simple simulation of 35 randomly oriented 
clouds in which half were ellipsoids with aspect ratios of 10:1. Although this simulation also 
reproduces the observed projected surface areas of the 15 large clouds fairly well, the volume 
filling factor in this case is ~19%. 
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Fig. 19. — All dynamical cloud morphologies are overlaid and colorcoded as in Figure [T71 
The upwind heliocentric direction of the velocity vector for each cloud is indicated by the 
symbol, while the downwind heliocentric direction is indicated by the © symbol. The Sun 
is moving approximately antiparallel to the LISM clouds. The star symbols indicate sight 
lines of radio scintillation sources, and the series of three small clouds centered at / = 222° 
and b 



44° are the H I contours from iHeiles fc Trolandl (120031) , of the cold cloud recently 
identified to be within the Local Bubble by Meyer et al. ( 20061 ). 




Fig. 20. — Two subsets of dynamical cloud morphologies are overlaid and colorcoded as in 
Figure [T71 The upwind heliocentric direction of the velocity vectors are indicated by the ® 
symbols, while the downwind heliocentric direction is indicated by the © symbols. The top 
plot shows the projected morphological similarities shared by the LIC, G, and Mic clouds, 
indicating that clouds like the Mic could result from collisions of other clouds, in this case 
the LIC and G clouds. The bottom plot show s the clouds in close ang ular proximity to the 
cold cloud (shown here are H I contours from lHeiles fc Trolandl (l2003l ). which are the series 
of three small clump s centered at / = 222° and b = 44°), which was identified to be within 



the Local Bubble by lMeyer et al.l (120061 ). Note the alignment of the cold cloud matches well 
with the alignment of the high- velocity Gem Cloud in the same location. The compressional 
macroscopic motions between the surrounding warm dynamical clouds (e.g., the Gem Cloud 
with the slower moving Leo, Aur, and LIC clouds), shown here may be the origin mechanism 
of the observed cold material. 
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Fig. 21. — Distribution of the weighted mean values of turbulent velocity and depletion 
of iron {D{Fe) and magnesium {D{Mg)) for all clouds with more than one sight line having a 
turbulent velocity measurement. The errors are the dispersion about the weighted mean. The 
dashed line is a weighted minimum linear fit to the data. A clear correlation exists between 
cloud turbulent velocity and depletion for both elements, with a hnear correlation coefficient 
r = 0.69 and a probability that the distribution could be drawn from an uncorrelated parent 
population of only Pc = 1.7% for iron and r = 0.73 and = 1.2% for magnesium. It is 
likely that regions of high turbulence result from the dynamical interactions of clouds, which 
in turn produce shocks and heat any dust, returning metal ions from the dust to the gas 
phase. 
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Fig. 22. — The distribution of cloud velocity differences between all 15 LISM clouds are 
shown here for all sight lines that are predicted to traverse multiple clouds based on the 
spatial distribution of clouds shown in Figure [191 The predicted velocity components of 
all 15 clouds were calculated for a uniform sample of hypothetical lines of sight over the 
full sky. Velocity differences between the LIC and G clouds in directions in which both are 
observed are indicated by the red histograms. The mean thermal sound speed is ~8 km s~^. 
A significant fraction of possible velocity differences between LISM clouds include velocities 
greater than the thermal sound speed. Distributions in the radial (top), and transverse [l, 
middle; b, bottom) Galactic directions are shown. The blue histogram (scaled by a factor of 
50) represents the macroscopic velocity differences for LISM material near the Leo cold cloud. 
Significant compressional velocities of the warm LISM material in th e radial direc t ion rn ay 
be a mechanism for the origin of the cold dense material observed by iMeyer et al.l (120061 ) . 
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Fig. 23. — Comparison of LIC and G sight line physical properties as a function of angular 
distance from the LIC downwind direction. Sight lines through each cloud are distinguished 
by color (red: LIC; pale pink: G) and symbol (circle: LIC; square: G). Since the G Cloud is 
in the upstream direction, all G Cloud sight lines are at high A^. Except for one LIC value, 
all temperatures through the LIC are higher than the G Cloud temperatures. A correlation 
with angle is evident in the depletion of magnesium (r = 0.63, Pc = 0.012%) and iron 
(r = 0.80, Pc = 0.0041%). 
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Fig. 24. — Similar to Figure [23J but the comparison of LIC and G sight hne physical prop- 
erties is shown as a function of hydrogen column density. Correlations are seen between 
hydrogen column density and depletion of magnesium (r = —0.53, Pc = 0.12%) and iron 
(r = -0.65, Pc = 0.15%). 
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Table 1. LIC Cloud Sight Line Properties 



star 


HD 


d 


V 






logiV(H I) 


T 






0(Mg) 


Other 


Name 




(pc) 


(km s~l) 






(cm-^) 


(K) 


(km s~^) 






Glouds*^ 


a CMa A 


48915 


2.6 


18.55 


0, 


.12 


17.2 


800oi™go 


0.5 ± 0.3 


— U. 'o_o 21 


— u-o4_Q 22 




a CMa B 


48915B 


2.6 


17.60 


0, 


.51 


17.6 


8000+?°° 
— 1000 


0.5 ± 0.3 








e Eri 


22049 


3.2 


18.73 


0. 


.70 


17.8 




2 03+° " 
■'■"'^-0.45 


— 1.03 it 0.17 


—0.81 it 0.12 




a CMi 


61421 


3.5 


19.76 


0. 


.78 


17.9 


6710+^^° 


1.21+0-55 
-0.45 


— 1.13 ± 0.07 


— 1.31 it 0.10 


Aur 


e Ind 


^uyiuu 


3.6 


— 9.20 


0. 


.02 


18.0 


<9850 








Vel 


T Get 


10700 


3.6 


12.34 


0. 


.20 


18.0 


<6700 








G, Blue 


40 Eri A 


26965 


5.0 


21.73 


0. 


.68 


17.8 


8120 lb 450 


n E+1.2 
0-5Io.5 




—0.90 it 0.10 




-q Gas A 


4614 


6.0 


11.18 


0, 


.67 












a PsA 


216956 


7.7 


-5.87 





,42 












Mic 


a Lyr 


172167 


7.8 


-12.90 


0, 


,30 












G 


Ori 


39587 


8.7 


23.08 





,38 


17.8 


Tnnn~b730 
7000_ggo 


„ ,0+0.15 
'=-3»_o.l7 


— i.Uo it 0.U7 


—0.74 it o.uy 




5 Eri 


23249 


9.0 


19.60 


0, 


,07 


17.9 


<8900 








Blue, Hyade 


K.^ Get 


20630 


9.2 


20.84 


0, 


,31 


17.6 


520011?™ 


■'■°*-0.32 


-1.01 ± 0.17 


-0.68 it 0.16 




CF UMa 


103095 


9.2 


2.05 


1 


,19 














f3 Gem 


()2,)()9 


10.3 


19.6,) 




,60 


18.0 


'"""'-1.100 


1 (;7+''-27 

^■''^-0,32 


1 -t- n 1 


1 f 1 '1 4- '1/1 




EP Eri 


17925 


10.4 


19.50 


0, 


,88 


18.0 


8100 ± 1300 


2.46 ± 1.45 




—0.48 ± 0.42 


G, Blue 




142806 


11.1 


— 18.19 





.93 












IVIic Leo 






























12 2 


17 50 

■ 


1 


25 

■ 




<9900 








ur 


a Aur 


34029 


12.9 


21.48 





.82 


18.2 


,,_„„+ 1400 
67001 1300 


1 <,„+0.32 
^■°°-0.39 


— 1.21 it 0.06 


—0.95 it 0.10 




HR 8 


166 


13.7 


6.50 


0, 


.44 


18.3 


< 10700 








Eri, Hyades 


a Oph 


159561 


14.3 


— 22.57 


1 


.68 














72 Her 


157214 


14.4 


-15.51 


0, 


.33 














a Gep 


203280 


15.0 


0.20 


0, 


.64 














CT Boo 


128167 


15.5 


—2.68 





.63 














99 Her 


165908 


15.7 


— 17.43 


1 


.64 












G 


/3 Gas 


432 


16.7 


9.16 


0, 


.69 


18.2 


9760i|™ 


0+ 


-1.28 it 0.16 


-1.28 it 0.17 




DX Leo 


82443 


17.7 


11.00 


0, 


.22 


17.7 


8540 ± 850 


1.78 ± 1.77 




-0.83 it 0.20 




Eri 


23754 


17.9 


16.99 


0, 


,04 














V368 Gep 


220140 


19.7 


6.00 


0, 


.19 


18.0 


12050+^^° 


0+^ ° 




-0.78 it 0.17 




a Tri 


11443 


19.7 


17.89 


1 


.84 


18.1 


7700+?],™ 
' ' "" — 2600 


0.0+ii 

— 0.0 


— 1.16 it 0.23 


— 1.44 it 0.15 




HR 4345 


97334 


21.7 


4.30 





,14 


17.8 


<8700 








NGP 


PW And 


1405 


21.9 


8.50 


0, 


,13 


18.1 


11300+11™ 
— 1800 


0.0+1 = 
— 0.0 




-0.48 it 0.52 


Eri, Hyades 


•^AO 1*^6 111 


73350 


23.6 


12.00 





,,50 


18.2 


<11400 








G Aur 


O iVlcl 


1 nR^Ql 
xuuoyj- 


25 


3 80 




















40183 


25 2 


22 30 




74 














■7 UMa 




25.6 


4.40 




.64 














X And 


222107 


25.8 


6.50 





.03 


18.6 


<13100 








H!yades 


(7 Get 


15798 


25.8 


15.99 





,81 






































nrv. ouu 


17948 


26 


l"" 10 





31 


















97 l" 


9 fio 


















MT? 1 riQQ 




29 


21 90 




88 


■ 


/yuu it louu 


1.10 i \j.'±i 




— 1 . zo in u. lu 




fl Poo- 


210418 


29 6 


4 20 




30 












Eri 


OrtW 0OUIO 


157466 


29 8 


j^g Q2 


Q 


96 
















'if ''38' 


'30^ 


13 O"* 
















Hyades 


5 Gas 








1 


06 
















209952 


31.1 


— 10.93 





.58 












Vel 


CK Lac 


213558 


31.4 


3.50 


1 


.11 






, „^+0.18 






Hyades 


DK UMa 


82210 


32.4 


9.41 


0, 


.06 


17.9 


6750±240 


-1.10 ± 0.07 


-0.87 ± 0.09 




e Gru 


216789 


39.7 


-7.30 


0, 


.43 












Vel 


SAO 76593 


27808 


40.9 


23.10 


0, 


.17 












Aur 


SAO 93981 


28568 


41.2 


23.90 


0. 


.39 


18.0 


<11700 






-0.97 ± 0.18 


Aur 


SAO 111879 


28736 


43.2 


21.60 


1 


.51 












Aur 


101 Tau 


31845 


43.3 


22.40 


1 


,44 












Aur 


SAO 94033 


29225 


43.5 


22.50 


1 


,10 












Aur 


SAO 93982 


28608 


43.6 


23.20 


0, 


.20 












Aur 


SAO 76683 


29419 


44.2 


23.20 


0, 


.26 












Aur 


V993 Tau 


28205 


45.8 


23.30 


0, 


.14 


18.0 


<11200 






-0.94 lb 0.29 


Aur 


SAO 93963 


28406 


46.3 


22.10 


1, 


,37 












Aur 


SAO 76609 


28033 


46.4 


23.60 


0, 


.28 


18.2 


< 12900 






-0.79 ± 0.36 


Aur 


V471 Tau 




46.8 


20.90 


0, 


,52 


18.2 


<11000 








Aur, Hyades 


SAO 93945 


28237 


47.2 


22.40 





,98 












Aur 


SAO 93831 


26784 


47.4 


23.00 


0, 


.05 












Aur 


7 Aqr 


212061 


48.4 


-4.62 


0, 


.76 
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Table 1 — Continued 



star 


HD 


d V 


CT» log JV(H I) 


T 






D(Mg) 


Other 


Name 


# 


(pc) (km s~^) 


(cm-2) 


(K) 


(km 






Clouds'" 



SAO 56530 


21847 


48.9 


21.10 


0.39 








SAO 93885 


27561 


51.4 


22.20 


1.11 






Aur 


K And 


222439 


62.0 


7.60 


1.06 






Hyades 


S Cyg 


186882 


62.4 


-9.60 


1.47 






G 


SAO 93913 


27848 


53.4 


22.40 


0.98 






Aur 


HR 1608 


32008 


54.7 


21.60 


0.07 


17.8 




Blue 


45 Aur 


43905 


57.0 


18.47 


0.13 








C Peg 


214923 


64.1 


-2.20 


0.34 






Eri 


t Cap 


203387 


66.1 


-12.06 


0.06 


17.9 


12900l||™ l-SSl o.lg -0.85 ± 0.36 -0.22 ± 0.43 




77 Aur 


32630 


67.2 


23.00 


0.94 








G191-B2B 




68.8 


19.19 


0.09 


18.2 


6200j:j3™ 1.78+o l° -1.18 ±0.09 -0.97 ± 0.10 




L Oph 


152614 


71.7 


-21.60 


0.13 








Feige 24 




74.4 


17.60 


0.17 


18.1 




Hyades 


7 Ori 


36468 


74.6 


25.40 


1.88 






Aur 


»(T = duo - 




where we 


have imposed a minimum CTy of 1 k 


m s~^for all high resolution data, and a minimum (Ty of 3 km 


s~^for all 



medium resolution data. 

^ Other dynamical clouds for which the observed velocity is within 3cr of the predicted value and the sight line is closer than 10° to the 
cloud. These are less likely, but possible dynamical assignments for these sight lines. 



Table 2. G Cloud Sight Line Properties 



star 


HD 


d 


V 




log Ar(H I) 


T 


( 


i3(Fe) 


D(Mg) 


Other 


Name 


# 


(pc) 


(km 




(cm-2) 


(K) 


(km s-1) 






Clouds'" 


a Ccn B 


128621 


1.3 


-18.14 


0.02 


17.6 


5500l||5j 


-'■•^'-0.41 


-0.53 ± 0.14 


-0.38 ± 0.14 




a Ccn A 


128620 


1.3 


-18.45 


0.33 


17.6 


5100+1200 


, ,,+0.33 
'-■^'--0.49 


-0.61 ± 0.06 


-0.42 ± 0.11 




70 Oph 


165341 


5.1 


-26.50 


0.07 


17.8 


2700+3000 
■^'00-2300 




-0.38 ± 0.09 


0.22 ± 0.20 


Mic 


36 Oph A 


155886 


5.5 


-28.40 


0.17 


17.8 


5870 ± 560 


„ ,„+0.46 


-0.61 ± 0.26 


-0.29 ± 0.18 




a Lyr 


172167 


7.8 


-16.86 


1.33 










Mic 


AB Dor 


36705 


14.9 


5.19 


0.15 












Blue, Dor 


LQ Hya 


82558 


18.3 


6.50 


0.23 


18.8 


syOO+IS^OO 
o '00-5700 


2.52 ± 1.26 




-1.38 ± 0.31 


Lie, Aur 


a Hyi 


12311 


21.9 


4.90 


0.76 












Get 


5 Vel 


74956 


24.4 


1.30 


1.18 














C Aql 


177724 


25.5 


-21.30 


0.30 












Eri 


CT Get 


15798 


25.8 


21.85 


0.39 














Eri 


18978 


26.4 


20.90 


0.12 














HR 4023 


88955 


31.5 


-1.70 


0.36 














SAO 68491 


184499 


32.0 


-14.28 


0.41 












Lie 


/3 Gar 


80007 


34.1 


-4.32 


1.66 














fj Vel 


93497 


35.3 


-4.38 


0.41 


18.6 


< 10600 








Get 


Lib 


136742 


49.1 


-26.90 


0.61 














o Ser 


160613 


51.5 


-29.00 


0.33 












Oph 


SAO 159459 


140283 


57.3 


-28.55 


0.02 












Gem 


Scr 


156928 


59.3 


-27.70 


1.49 












Mic 


HD 141569 


141569 


99.0 


-28.70 


0.22 












Gem 


V = (1^0 - 


Vi,\)/a^, 


where we 


have imposed 


a minimum o"^ of 1 


km s~lfor all 


high resolution 


data, and a minimum cr^ of 3 km s ^ for all 



medium resolution data. 

Other dynamical clouds for which the observed velocity is within 3a of the predicted value and the sight line is closer than 10° to the cloud. 
These are less likely, but possible dynamical assignments for these sight lines. 
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Table 3. Blue Cloud Sight Line Properties 



Star 


HD 


d 


V 




logiV(H I) 


T 


e 


D(Fe) 


D(Mg) 


Other 


Name 


# 


(pc) 


(km s-^) 




(cm-2) 


(K) 


(km 






Clouds'" 


a CMa A 


48915 


2.6 


12.70 


0.04 


17.2 


•^"""-1000 


2.7 ± 0.3 


-0.95lSi5 


-0 77+° " 




a CMa B 


48915B 


2.6 


11.70 


0.63 






2.7 ± 0.3 








BP Bri 


17925 


10.4 


9.00 


0.65 












Lie 


C Dor 


33262 


11.7 


8.41 


0.25 


17.8 


--„„+2300 
' '^"-2100 


2 34+0.38 


-0.52 ± 0.30 


-0.05 ± 0.27 


G 


HR 2225 


43162 


16.7 


14.00 


0.44 


17.9 


< 10400 








Lie 


(3 Pic 


39060 


19.3 


10.12 


0.81 














HR 2882 


59967 


21.8 


8.40 


0.47 


18.6 


<15100 










CT Get 


15798 


25.8 


9.74 


0.38 














/3 Car 


80007 


34.1 


3.67 


0.24 












Vel, Cet 


a Car 


45348 


95.9 


8.19 


0.68 












Dor 



^(T — (1^0 ~ 'V-k\)/^v> where we have imposed a minimum of 1 km s "'^for all high resolution data, and a minimum (T^, of 3 km s "'^for 
all medium resolution data. 

^Other dynamical clouds for which the observed velocity is within 3(7 of the predicted value and the sight line is closer than 10° to the 
cloud. These are less likely, but possible dynamical assignments for these sight lines. 



Table 4. Aql Cloud Sight Line Properties 



star 


HD 


d 


V 


ct'' logJV(Hl) T J D(Fe) 




Other 


Name 


# 


(pc) 


(km s-') 


(cm-2) (K) (km s-') 




Clouds'" 


61 Cyg A 


201091 


3.-5 


-:-i.()0 


1.28 17.8 6850 ± 880 2.08 ± 0.64 


-0.95 ± 0.22 


Lie 


a Aql 


187642 


5.1 


-20.90 


0.14 17.9 12600 ± 2400 0.631°;^° -0.96 ± 0.54 


-0.69 ± 0.55 


G, Bri 


70 Oph 


165341 


5.1 


-43.34 


0.92 17.1 3300 ± 2100 2.31 ± 0.37 


-0.67 ±0.16 




HR 6748 


165185 


17.4 


-29.20 


0.36 18.1 <15700 




G 


C Aql 


177724 


26.6 


-30.20 


0.79 




Oph 


\ Aql 


177756 


38.4 


-30.70 


1.19 




Mic 


BO Mic 


197890 


AAA 


0.00 


0.31 






i Cyg 


186882 


52.4 


-18.80 


0.71 




Mic 


a Del 


196867 


73.8 


-10.40 


1.01 






*CT = (\vo 


-•u*l)/<^„, 


where we have imposed a minimum (j^, of 1 km s "'for all high resolution data, and a minimum (Ty of 3 km 


s~'for all 


medium resolution data. 










''Other dynamical clouds for which the observed velocity is within 3(j of the predicted value and the sight line is 


closer than 10° to the cloud. 


These are less likely, but possible dynamical assignments for these sight lines. 












Table 5. Eri Cloud Sight Line Properties 






Star 


HD 


d 


V 


CT* logAr(Hl) T J £>(Fe) 


i3(Mg) 


Other 


Name 


# 


(pc) 


(km s-') 


(cm-2) (K) (km s-') 




Clouds'" 



61 Cyg A 


201091 


3.5 


-9.00 


0.27 


17.8 


6860 ± 880 


2.08 ± 0.64 




-0.36 ± 0.22 




a Aql 


187642 


5.1 


-17.10 


1.04 


17.9 


12300+2°°° 




-0.86 ± 0.27 


-0.67 ± 0.28 


G 


A Aql 


177756 


38.4 


-21.90 


0.35 












G 


a Peg 


218045 


42.8 


-0.90 


0.12 












Lie 


V376 Peg 


209458 


47.1 


-6.60 


0.00 


18.4 


< 15900 










V Peg 


220657 


53.1 


1.73 


0.13 


17.9 


1700+"°° 


3.93 ± 0.22 


-0.35 ± 0.08 


0.01 ± 0.12 


Lie 


a Del 


196867 


73.8 


-14.00 


0.14 












e Aql 


191692 


88.0 


-17.70 


0.15 












Aql 



= {\vq — |}/cTu! where we have imposed a minimum cr^ of 1 km s ^for all high resolution data, and a minimum o"^ of 3 km s ""^for all 
medium resolution data. 

^Other dynamical clouds for which the observed velocity is within 3ct of the predicted value and the sight line is closer than 10° to the 
cloud. These are less likely, but possible dynamical assignments for these sight lines. 



-58- 



Table 6. Aur Cloud Sight Line Properties 



star 


HD 


d 


V 




log JV(H I) 


T 






D(Fe) D(Mg) 


Other 


Name 


# 


(pc) 


(km s-^) 




(cm-2) 


(K) 


(km s~ 






Clouds'' 


a CMi 


61421 


3.5 


23.00 


1.03 


17.6 




i.2iiS: 


35 
.45 


-1.13 ±0.07 -0.79 ±0.10 




LQ Hya 


82558 


18.3 


14.00 


0.19 












Lie, G 


a Tau 


29139 


20.0 


20.62 


1.29 


17.8 










Lie 


SAO 93801 


26345 


43.1 


21.10 


0.51 












Lie 


SAO 93973 


28483 


50.2 


20.80 


0.31 












Lie 


7 Crv 


106625 


50.6 


-2.00 


0.32 














45 Aur 


43905 


57.0 


11.87 


0.60 














?7 Aur 


32630 


67.2 


17.44 


0.87 














HR 2324 


45320 


70.1 


26.00 


1.18 















^(7 — {\vq — Tii)/cr,;, where we have imposed a minimum of 1 km s "'^for all high resolution data, and a minimum of 3 km s "'^for 
all medium resolution data. 

^ Other dynamical clouds for which the observed velocity is within 3tT of the predicted value and the sight line is closer than 10° to the 
cloud. These are less likely, but possible dynamical assignments for these sight lines. 



Table 7. Hyades Cloud Sight Line Properties 



star 


HD 


d 


V 






logiV(H I) 


T 


f 


D(Fe) 


D(Mg) 


Other 


Name 


# 


(pc) 


(km s-^) 






(cm-2) 


(K) 


(km 






eiouds'' 


BV Lac 




5.0 


7.30 


0. 


,12 


18.0 


<8300 








Lie, Eri 


Get 


20630 


9.2 


13.36 


1. 


11 


17.6 


3600 + ??™ 


, , 7+0.34 


-1.29 ± 0.26 


-0.79 ± 0.19 




a Tri 


11443 


19.7 


13.66 


0. 


16 


17.8 


89«<)li™ 


1 

-^"'-1.3 


-0.97 ± 0.23 


-0.58 ± 0.15 


LIG 


HR 1099 


22468 


29.0 


14.80 


0, 


,63 


17.6 


8800l»?°„ 




< -0.82 


-1.36 ± 0.0.09 




SAO 93801 


26345 


43.1 


13.60 


0. 


08 














SAO 111879 


28736 


43.2 


13.40 


0, 


19 














SAO 76683 


29419 


44.2 


12.30 


0, 


,54 














V993 Tau 


28206 


45.8 


14.80 


0. 


,84 














SAO 93831 


26784 


47.4 


15.50 


1. 


,62 














SAO 93886 


27661 


51.4 


14.40 


0. 


,86 














V Peg 


220667 


53.1 


8.80 


0. 


,51 


17.4 


1000+^^°° 


^ -0.63" 


-0.87 ± 0.45 


-1.07 ± 0.30 




45 Aur 


43906 


57.0 


8.14 


0. 


,84 














ri Aur 


32630 


67.2 


10.70 


0. 


,02 














G191-B2B 




68.8 


8.61 


0. 


,74 


17.4 


4400+^*°° 


o ,7+0.37 
' -0.39 


0.05 ± 0.11 


0.87 ± 0.54 





^o" = (^\vq — v^\) / <y^ , where we have imposed a minimum of 1 km s "'^for all high resolution data, and a minimum <t^ of 3 km s "'^for all 
medium resolution data. 

'^Other dynamical clouds for which the observed velocity is within 3(T of the predicted value and the sight line is closer than 10° to the cloud. 
These are less likely, but possible dynamical assignments for these sight lines. 
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Tablc 8. Mic Cloud Sight Line Properties 



star 


HD 


d 


V 






log JV(H 1) 


T 


€ 


i3(Fe) 


D(Mg) 


Other 


Name 


# 


(pc) 


(km s~^) 






{cm-2) 


(K) 


(km s~l) 






Clouds'" 


a Aql 


187642 


5.1 


-25.02 


0, 


,22 


17.5 


12B00i=™g 


1 4+°-'' 


-0.61 ± 0.23 


-0.38 ± 0.24 




a PsA 


216956 


7.7 


-10.64 





()7 














CK Lyr 


172167 


7.8 


-19.40 


0, 


,19 














AU Mic 


197481 


9.9 


-21.45 


0, 


,05 


18.2 


8700 ± 1200 


4.30 ± 0.93 




-0.55 ± 0.19 




a Oph 


159561 


14.3 


-26.23 


0, 


,05 












G 


99 Her 


165908 


15.7 


-22.90 


0, 


.60 














a CrB 


139006 


22.9 


-17.40 


0, 


,21 












Oph, Leo 


/3 Get 


4128 


29.4 


1.63 


0, 


,15 


16.9 


12400 ± 2800 


2.29 ± 0.44 




0.23 ± 0.09 




SAO 68491 


184499 


32.0 


-19.73 





,74 














A Aql 


177756 


38.4 


-26.50 


1 


,38 












G 


/3 Sor 


141003 


46.9 


-20.70 


0, 


,47 












Lie, Oph, Leo 


A Oph 


148857 


50.9 


-24.80 


0, 


,07 












Lie, Oph 


5 Cyg 


186882 


52.4 


-16.30 


0. 


.27 














I Cap 


203387 


66.1 


-20.48 


0, 


.41 


18.1 


11700+*^°° 


o 09+0.37 
3.»^_0.44 


-1.22 ± 0.23 


-0.47 ± 0.21 




e CrB 


138749 


95.3 


-16.70 


0, 


.29 












Oph 



= (|vo — "^isY) / f^v 1 where we have imposed a minimum <Ty of 1 km s ^for all high resolution data, and a minimum Uy of 3 km s ^for all medium 
resolution data. 

'^Other dynamical clouds for which the observed velocity is within 3(t of the predicted value and the sight line is closer than 10° to the cloud. These 
are less likely, but possible dynamical assignments for these sight lines. 



Table 9. Oph Cloud Sight Line Properties 



star 


HD 


d 


V 




log ]V(H I) 


T 




■D(Fe) 


£>(Mg) 


Other 


Name 


# 


(pc) 


(km s-"-) 




{cm-2) 


(K) 


(km g-"-) 






Clouds'' 


70 Oph 


166341 


5.1 


-32.53 


1.30 


17.5 


1700+^^™ 
^'""-1700 


3.3 ± 1.1 




-0.84 ± 0.34 




7 Ser 


142806 


11.1 


-22.12 


0.28 










Mic 


ct Oph 


159561 


14.3 


-28.40 


1.17 












G, Mic 


72 Her 


157214 


14.4 


-25.95 


0.32 












NGP 


7 Oph 


161868 


29.1 


-29.90 


0.62 












G, Mic 


SAO 68491 


184499 


32.0 


-27.85 


0.02 















^cr — {\vq — v^\)/a^). where we have imposed a minimum of 1 km s "'^for all high resolution data, and a minimum cr^, of 3 

km s~"'"for all medium resolution data. 

'^^ Other dynamical clouds for which the observed velocity is within 3cr of the predicted value and the sight line is closer than 10° 
to the cloud. These are less likely, but possible dynamical assignments for these sight lines. 
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Table 10. Gem Cloud Sight Line Properties 



star 


HD 


d 


V 




log iV(H I) 


T 




D(Fe) D(Mg) 


Other 


Name 


# 


(PC) 


(km 




(cm-2) 


(K) 


(km 




Clouds'" 


{ Boo A 


131156 


6.7 


-17.69 


0.33 


17.9 


5310 ± 830 


1.68 ± 0.23 


-0.92 ± 0.10 


Mic, NGP 


p Gem 


62509 


10.3 


31.84 


0.00 


17.8 


fiinn + 3100 


1 93+0-^9 


-1.29 ±0.12 -1.19 ±0.14 




ct Boo 


124897 


11.3 


-13.89 


0.50 










NGP 


<S Crv 


108767 


26.9 


-0.50 


0.76 










Aur 


HR 4803 


109799 


34.6 


-0.37 


1.02 










Aur 


a Gem 


62044 


37.5 


32.26 


0.55 


17.7 


_2„„+iooo 

'^UU_1200 


0+1-^ 


-1.18 ± 0.14 




Ser 


141003 


46.9 


-23.30 


1.61 










Oph, NGP 


g Lib 


135742 


49.1 


-23.60 


0.17 












7 Crv 


106625 


50.6 


1.60 


0.73 












Cen 


129685 


63.5 


-16.90 


0.36 













^CT = (1^0 — ''^■tt \ ) / <^v > where we have imposed a minimum ax, of 1 km s "'^for all high resolution data, and a minimum of 3 km s "'^for all 
medium resolution data. 

^Other dynamical clouds for which the observed velocity is within 3(7 of the predicted value and the sight line is closer than 10° to the cloud. 
These are less likely, but possible dynamical assignments for these sight lines. 



Table 11. NGP Cloud Sight Line Properties 



star 


HD 


d 


V 




log JV(H I) 


T 




C(Fe) 


0(Mg) 


Other 


Name 


# 


(pc) 


(km 




(cm-2) 


(K) 


(km 






Clouds'" 


61 Vir 


115617 


8.5 


-16.50 


0.01 


17.9 


<8500 








Leo 


a Oph 


159561 


14.3 


-32.74 


0.84 














X Her 


142373 


15.9 


-12.90 


0.76 


18.2 


< 10200 








Mic 


c Cen 


115892 


18.0 


-18.20 


0.31 












G 


SAO 28753 


116956 


21.9 


3.20 


0.57 


18.2 


< 12600 








Lie 


I Leo 


99028 


24.2 


4.99 


0.65 












Leo 


77 UMa 


120315 


30.9 


-3.02 


0.11 


17.9 


8900_23oj, 




-0.93 ± 0.06 


-0.78 ± 0.10 




L Dra 


137759 


31.3 


-6.81 


1.74 












Lie 


HZ 43 




32.0 


-6.52 


0.40 


18.0 


7500+^"^°° 


1 7+0.8 
^■'-1.7 


-1.42 ± 0.11 


-1.08 ± 0.17 




CVn 


112413 


33.8 


-1.90 


0.77 












Lie 


HR 4803 


109799 


34.6 


-11.41 


0.55 














Cen 


129685 


63.5 


-26.30 


0.33 












G 


GD 153 




70.5 


-5.04 


0.44 


17.9 


7000^^2^ 


1 2+1-^ 


-1.18 ± 0.13 






31 Com 


111812 


94.2 


-3.37 


0.03 


18.0 


S200+_\°Z 


0+1° 


-0.85 ± 0.15 


-0.98 ± 0.14 




CrB 


138749 


95.3 


-19.90 


0.72 












Oph, Gem 



^(T = (\vq — u^|)/(Tu, where we have imposed a minimum a^j of 1 km s "'^for all high resolution data, and a minimum o"^ of 3 km s "'^for all 
medium resolution data. 

Other dynamical clouds for which the observed velocity is within 3ct of the predicted value and the sight line is closer than 10° to the cloud. 
These are less likely, but possible dynamical assignments for these sight lines. 
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Table 12. Leo Cloud Sight Line Properties 



Star 


HD 


d 


V 




log Af(H I) 


T 5 D(Fe) D(Mg) 


Other 


Name 


# 


(pc) 


(km s-^) 




(cm-2) 


(K) (kms-1) 


Clouds'" 


/3 Leo 


102647 


11.1 


0.51 


0.05 






NGP 


HR 4657 


106616 


22.6 


-3.40 


0.42 


18.6 


<27500 


Aur, Gem, NGP 


a Leo 


87901 


23.8 


10.50 


0.76 






Lie 


i Leo 


99028 


24.2 


1.97 


0.50 








HR 4803 


109799 


34.6 


-7.52 


0.40 








Ser 


141003 


46.9 


-16.70 


1.05 






Lie 


Oi Vir 


116658 


80.4 


-11.99 


0.80 








"ct = (|l>o - 


- 


, , where 


we have imp 


osed a 


minimum of 


1 km s^^for all high resolution data, and 


a minimum of 3 


km ^ for all 


medium 


resoluti 


on data. 










^ Other dynamical clouds for which the observed velocity is withi 


n 3(T of the predicted value and the sight li 


ne is closer than 10° 



to the cloud. These are less likely, but possible dynamical assignments for these sight lines. 



Table 13. Dor Cloud Sight Line Properties 



Star 
Name 


HD 
# 


d 
(pc) 


V 

(km s-1) 




log JV(H I) 
(cm-2) 


T 
(K) 


(km s-1) 


-D(Fe) 


D(Mg) 


Other 
eiouds'' 


C Dor 
t"' Eri 
e Gru 
a Eri 


33262 
23754 
215789 
10144 


11.7 
17.9 
39.7 
44.1 


13.90 
41.40 
12.20 
21.20 


0.35 
0.12 
0.66 
0.24 


18.1 


7000+2^"'' 

/■UUU_3QQQ 


5 47+0-39 


-0.80 ± 0.30 


-0.66 ± 0.31 





^(j — (|tJQ — iJ*!)/ (7^ , where we have imposed a minimum ay of 1 km s "'^for all high resolution data, and a minimum Cy of 3 km s "'^for 
all medium resolution data. 

Other dynamical clouds for which the observed velocity is within 3<7 of the predicted value and the sight line is closer than 10° to 
the cloud. These are less likely, but possible dynamical assignments for these sight lines. 



Table 14. Vel Cloud Sight Line Properties 



star 


HD 


d 


V 


a* 


logJV(H I) 


T 


« -D(Fe) 


i3(Mg) 


Other 


Name 


# 


(pc) 


(km s-1) 




(cm-2) 


(K) 


(km s-1) 




eiouds'' 


AB Dor 


36706 


14.9 


14.22 


0.11 










Get 


SAO 254993 


254993 


20.5 


-13.10 


1.24 


18.8 


< 23000 






Lie, G 


a Hyi 


12311 


21.9 


9.80 


0.31 












.5 Vel 


74956 


24.4 


11.80 


0.18 












a Eri 


10144 


44.1 


11.00 


0.45 












BO Mic 


197890 


44.4 


-24.20 


0.28 


18.3 


10600 ± 2700 


3.48 ± 1.87 


-0.03 ± 0.72 


Mic 


a Pav 


193924 


66.2 


-19.60 


0.63 












"■a = (\vo - 


V* |)/cr^, whe 


re we 


have imposed 


a minimum (7^ of 1 km 


s "'^for all high resolution data, and a 


minimum o"^ of 3 k 


m s~^for 



all medium resolution data. 



^ Other dynamical clouds for which the observed velocity is within 3a of the predicted value and the sight line is closer than 10° to the 
cloud. These are less likely, but possible dynamical assignments for these sight lines. 
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Table 15. Get Cloud Sight Line Properties 



star 
Name 


HD 

# 


d 
(pc) 


V 

(km 




logJV(H I) 
(cm-2) 


T 
(K) 


S D(Fe) 

(km 


C(Mg) 


Other 
Clouds'' 


cr Boo 


128167 


15.5 


-12.90 


1.10 












6 Vel 


74956 


24.4 


15.60 


0.79 












/3 Get 


4128 


29.4 


9.14 


0.53 


18.5 


6300 ± 2900 


1.31 ± 0.76 


0.21 ± 0.33 


Lie 


IIR 402.'-! 


88')r,r, 


■■VLSt 


10. rM 


0,48 












(> Eli 


1(1144 


44,1 


7.(i0 


1,80 










C. Vol 



*CT = (|vo — v*|)/(Ti,, where we have imposed a minimum a^, of 1 km s ■'^for all high resolution data, and a minimum cr^ of 3 
km s~"'^for all medium resolution data. 

^Other dynamical clouds for which the observed velocity is within 3a" of the predicted value and the sight line is closer than 10° 
to the cloud. These are less likely, but possible dynamical assignments for these sight lines. 
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Table 16. LISM Cloud Heliocentric Velocity Vectors 



Cloud 


Number of 


Vo 


^0 


&o 


9 

xl 


Name 


Sight Lines 


/I 1 \ 
(km s j 


(°) 


(°) 




Lie 


7Q 


OQ C4 _|_ n on 




„i Q c: _L 

±0.0 —1— 0.0 


2.2 


G 


21 


29.6 ± 1.1 


184.5 ± 1.9 


-20.6 ± 3.6 


1.3 


Blue 


10 


89 ± 89 


205 5 ± 4 3 


-21.7 ± 8.3 


2.4 


Aql 


9 




1 07 A 1 1 C 


_cn 1 1 A 

UL/.O —1— J-.v7 


2.6 


Eri 


Q 

o 


24.1 ± 1.2 


1 ofi 7 + 2 1 


_i 7 7 _|_ 2 fi 


3 


A.ur 


9 


25 22 + n 81 


21 2 n + 2 4 


-16.4 ± 3.6 


2.1 




14 


1 4 CO + n Qi 

J.4:.U-7 —1— VJ.OJ. 


1 «4 2 + Q 4 


-42.8 ± 6.1 


1 3 


Mic 


15 


28.45 ± 0.95 


203.0 ± 3.4 


-03.3 ± 2.3 


0.5 


Oph 


6 


32.25 ± 0.49 


217.7 ± 3.1 


+00.8 ± 1.8 


3.9 


Gem 


10 


36.3 ± 1.1 


207.2 ± 1.6 


-01.2 ± 1.3 


1.7 


NGP 


15 


37.0 ± 1.4 


189.8 ± 1.7 


-05.4 ± 1.1 


3.8 


Leo 


7 


23.5 ± 1.6 


191.3 ± 2.8 


-08.9 ± 1.8 


1.5 


Dor 


4 


52.94 ± 0.88 


157.3 ± 1.5 


-47.93 ± 0.63 


0.8 


Vel 


7 


45.2 ± 1.8 


195.4 ± 1.1 


-19.1 ± 1.0 


0.8 


Get 


5 


60.0 ± 2.0 


197.11 ± 0.56 


-08.72 ± 0.50 


8.9 


LIG'' 


9 


25.7 ± 0.5 


186.1 


-16.4 




LIG'^ 


16 


26 ± 1 


186 ± 3 


-16 ± 3 




LIG'^ 


63 


24.20 ± 1.05 


187.0 ± 3.1 


-13.5 ± 3.0 


2.1 


G^ 




29.4 


185.5 


-20.5 




Helio'^ 




26.24 ± 0.45 


183.4 ± 0.4 


-15.9 ± 0.4 




(LIG+G)/2° 




26.74 ± 0.71 


185.7 ± 3.4 


-16.95 ± 3.6 





Lallement fc BertiiJ |l992h 
' iLallement et al.l |l995h 



'^LIC flow vector deleting the 16 lines of sight near the decelerated leading edge of 
the Lie in the direction of the Hyades Gloud. 

''Flow vector for interstella r helium gas in the hcliosphere. Temperature is 6303 
± 390 K: iMobius et all (|2004r ). See temperatures for individual dynamical clouds in 
Table m 

° Average of the LIG and G vectors. Average temperature of the LIG and G Glouds 
is 6500 ± 680 K. The in situ "Helio" measurement is closer to the average LIG and 
G temperature, than either cloud individually, see Table fT8l 
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Table 17. Unassigned Sight Line Properties 



star 


HD 


d 


V 


log JV{H I) 


T 




D(¥c) 


D{Mg) 


Possible 


Name 


# 


(pc) 


(km s-^) 


(cm-2) 


(K) 


(km s-^) 






Clouds* 


AD Leo 




4.7 


13.13 


18.5 


< 12300 








LIG, Leo 


77 Cas A 


4614 


6.0 


5.80 














Get 


20630 


9.2 


7.36 


17.4 


5800 ± 2700 


1.48 ± 0.92 




-1.61 ± 0.23 




/3 Leo 


102647 


11.1 


11.80 














AB Dor 


36705 


IAS 


19.38 














Bri 


23754 


17.9 


27.93 














Eri 


23754 


17.9 


8.82 














PW And 


1405 


21.9 


2.00 












Lie, Eri 


SAO 158720 


128987 


23.6 


-22.00 


18.1 


< 12400 








G, Gem, NGP, Leo 


.5 Her 


166164 


24.1 


-19.60 












LIG 


SAO 32862 


198084 


27.1 


-12.88 














HR, 1099 


22468 


29.0 


8.20 


17.2 


7100 ± 1400 


2.30 ± 0.26 




-1.21 ± 0.12 




7 Oph 


161868 


29.1 


-33.00 












Oph, NGP 


a And 


358 


29.8 


13.00 












Hyades 


ri UMa 


120315 


30.9 


2.60 


16.6 


r,+4400 


e e+0.9 

5-6Ii.i 




—0.78 ± 0.16 


NGP 


a Gru 


209952 


31.1 


-6.40 












Lie, Vel 


OL Gru 


209952 


31.1 


-21.40 












Lie, Vel 


o" Gem 


62044 


37.5 


21.77 


17.9 


8600 ± 1600 


2.46 ± 0.45 




-0.92 ± 0.14 


Lie 


HR 2298 


44769 


39.4 


7.00 














e Gru 


216789 


39.7 


6.80 














e Gru 


215789 


39.7 


-1.70 














e Gru 


215789 


39.7 


-12.30 














e Gru 


215789 


39.7 


-22.40 














SAO 93981 


28568 


41.2 


16.50 












Hyades 


SAO 111879 


28736 


43.2 


-4.30 














SAO 93982 


28608 


43.6 


15.80 












Hyades 


a Eri 


10144 


44.1 


18.90 












Dor 




28237 


47 2 


15 60 














P Lib 


135742 


49.1 


-33.70 














SAO 94162 


30738 


51.8 


20.30 












Aur 


K And 


222439 


52.0 


0.80 














Peg 


220657 


53.1 


-7.48 


17.4 


3600+4400 

JDUU_3(j(j(j 


1 7+0.7 
^■'-1.5 


-0.83 ± 0.2 


-0.54 ± 0.21 




SAO 93913 


27848 


53.4 


16.40 












Aur, Hyades 


HR 1608 


32008 


54.7 


6.10 


17.7 


<11500 








Blue, Hyades 


a Pav 


193924 


56.2 


-18.60 












Vel 


ri Aqr 


213998 


56.3 


-2.10 












Lie 


SAO 159459 


140283 


57.3 


-22.48 












Gem, Leo 


I Cap 


203387 


66.1 


-2.22 


18.4 


5500+1°^°° 
0000-5500 


3 7+°-'' 
-1.1 


-0.63 ± 0.37 


-0.28 ± 0.38 




77 Aur 


32630 


67.2 


6.80 














L Oph 


152614 


71.7 


-11.00 














t Oph 


152614 


71.7 


-30.80 












NGP 


Feige 24 




74.4 


3.10 


18.2 












7 Ori 


35468 


74.5 


16.00 














7 Ori 


35468 


74.5 


20.20 














7 Ori 


35468 


74.5 


27.80 














a Vir 


116658 


80.4 


-5.40 












Gem 


a Vir 


116658 


80.4 


-3.68 














e Aql 


191692 


88.0 


-22.60 












Mic 


r Her 


147394 


96.4 


-33.60 














r Her 


147394 


96.4 


-38.00 














HD 141569 


141669 


99.0 


-13.70 















^Although they do not fit the criteria for membership in any of the dynamical clouds, we list here possible membership where the observed velocity 
is within 3(t of a predicted cloud value and the sight line is closer than 10° to the cloud. 



Table 18. Summary of Cloud Properties 



Cloud 


# of 


Gentral Goord. 


Closest 


Surface Area 




# of 


(0 


# of 


(D(Fe)) 


# of 


(C(Mg)) 


# of 


Morphology 


Name 


Sight Lines 




6(°) 


Star (pc) 


(sq. deg) 


(K) 


Stars 


(km s-1) 


Sight Lines 




Sight Lines 




Sight Lines 




Lie 


79 


170 


-10 


2.6 


18270 


7500 ± 1300 


19 


1.62 ± 0.75 


19 


-1.12 ± 0.10 


12 


-0.97 ± 0.23 


21 


compact 


G 


21 


315 


+00 


1.3 


8230 


5500 ± 400 


5 


2.2 ± 1.1 


5 


-0.54 ± 0.11 


4 


-0.36 ± 0.35 


5 


compact 


Blue 


10 


250 


-30 


2.6 


2310 


3900 ± 2300 


3 


2.64 ± 0.16 


3 


-0.84 ± 0.27 


2 


-0.51 ± 0.49 


2 


compact 


Aql 


9 


40 


-05 


3.5 


2960 


7000 ± 2800 


3 


2.07 ± 0.64 


3 


(-0.96)" 


1 


-0.69 ± 0.21 


3 


compact 


Eri 


8 


70 


-20 


3.5 


1970 


5300 ± 4000 


3 


3.6 ± 1.0 


3 


-0.39 ± 0.19 


2 


-0.15 ± 0.30 


3 


compact 


Aur 


9 


210 


+ 10 


3.5 


1640 


(6710)" 


1 


(1.2)" 


1 


(-1.13)" 


1 


(-0.79)" 


1 


filamentary 


Hyades 


14 


180 


-20 


5.0 


1810 


6200 ± 3800 


6 


2.7 ± 1.2 


6 


-0.32 ± 0.62 


4 


-1.06 ± 0.47 


5 


filamentary 


Mic 


15 


40 


+ 15 


5.1 


3550 


9900 ± 2000 


4 


3.1 ± 1.0 


4 


-0.92 ± 0.43 


2 


-0.03 ± 0.40 


4 


filamentary 


Oph 


6 


45 


+25 


5.1 


1360 


(1700)" 


1 


(3.3)" 


1 







(-0.84)" 


1 


compact 


Gem 


10 


300 


+40 


6.7 


3300 


6000 ± 1100 


3 


1.63 ± 0.41 


3 


(-1.29)" 


1 


-1.05 ± 0.16 


3 


filamentary 


NGP 


15 


5 


+75 


8.5 


4020 


8000 ± 600 


4 


1.23 ± 0.43 


4 


-1.04 ± 0.23 


4 


-0.89 ± 0.15 


3 


compact 


Leo 


7 


270 


+55 


11.1 


2400 






















compact 


Dor 


4 


270 


-50 


11.7 


1550 


(7000)" 


1 


(5.5)" 


1 


(-0.80)" 


1 


(-0.65)" 


1 


compact 


Vel 


7 


300 


-45 


14.9 


2190 


(10600)" 
(6300)" 


1 


(3.5)" 
(1.3)" 


1 







(-0.03)" 
(0.21)" 


1 


compact 


Get 


5 


290 


-40 


15.6 


2270 


1 


1 







1 


filamentary 



^Those clouds with only one sight line with a physical measurement are indicated in parentheses. Since a weighted average is not possible, the listed properties should be considered uncertain. 
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